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Abstract 

This document is an extended version of the results presented in S. Dobberschiitz: Effec¬ 
tive Behaviour of a Free Fluid in Contact with a Flow in a Curved Porous Medium, SIAM 
Journal on Applied Mathematics, 2015. 

The appropriate boundary condition between an unconfined incompressible viscous fiuid 
and a porous medium is given by the law of Beavers and Joseph. The latter has been justified 
both experimentally and mathematically, using the method of periodic homogenisation. 
However, all results so far deal only with the case of a planar boundary. In this work, we 
consider the case of a curved, macroscopically periodic boundary. By using a coordinate 
transformation, we obtain a description of the flow in a domain with a planar boundary, 
for which we derive the effective behaviour: The effective velocity is continuous in normal 
direction. Tangential to the interface, a slip occurs. Additionally, a pressure jump occurs. 

The magnitude of the slip velocity as well as the jump in pressure can be determined with 
the help of a generalised boundary layer function. The results indicate the validity of a 
generalised law of Beavers and Joseph, where the geometry of the interface has an influence 
on the slip and jump constants. 
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1 Introduction 

A now classical result in the theory of homogenization states that, starting with the Stokes or 
Navier-Stokes equation, the effective fluid flow in a porous medium is given by Darcy’s law (see 
the works of Tartar in |SP80| , Allaire in |Hor97| and Mikelic [MikQlj ). When dealing with porous 
bodies inside another fluid, the boundary condition coupling the free fluid flow and the Darcy 
flow at the porous-liquid interface is of great interest. However, due to the different nature of 
the governing equations, the derivation of a ‘natural’ boundary condition is difficult: While the 
equation for the Darcy velocity consists of a second order equation for the pressure and a first 
order equation for the velocity, the system of equations governing the free fluid velocity (e.g. the 
Stokes or Navier-Stokes equation) is of second order for the velocity and of first order for the 
pressure. 

For an incompressible fluid, the flow in the direction normal to the interface has to be contin¬ 
uous due to mass conservation. However, additional conditions at the interface are not clearly 
available. 

From a mechanical point of view, Beavers and Joseph |BJ67| concluded by practical experi¬ 
ments that a .jump in the effective velocity appears in tangential direction. This jump is given 
by 

aK ^ {vp — vd) ■ T = {\7 ■ T, (1-1) 

where vp denotes the velocity of the free fluid, vp denotes the effective Darcy velocity in the 
porous medium and K is the permeability of the porous medium. The factor a is the so-called 
slip coefficient which has to be determined experimentally. Moreover, u and r are the unit normal 
and tangential vector with respect to the interface separating the porous medium and the free 
fluid. The Darcy velocity in the absence of outer forces for given fluid viscosity /i is given by 

vd = - -A Vp, 

where p denotes the pressure. Note that the condition mentioned above gives a relation between 
the velocity of the free fluid at the interface and the effective velocity inside the porous medium 
- it does not impose a condition on the actual fluid velocity inside the porous medium at the 
interface. Figure contains a schematic illustration. 
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Later, SafFman used a statistical model to derive the boundary condition of Beavers and 
Joseph. In |Saf71] . he argued that • r is of lesser order than the other terms and arrived at a 
jump given by 

vf ■ T = —K^ (V vpi') ■ T + (1-2) 

a 

Other boundary conditions were proposed as well: Ochoa-Tapia and Whitaker for example used 
the REV-method to obtain that the velocity and pressure as well as the normal stress are con¬ 
tinuous over the porous-liquid interface, but a jump appears in the tangential stress in the form 

{y{vD)v - V(u_f)i/) • T = PK~^Vd) ■ T. 

Here (vp) denotes the averaged free fluid velocity, which is given by a Stokes equation, and 
(vd) is the averaged velocity in the porous medium, which in this case fulfills a Darcy law with 
Brinkman correction, 

(vd) = - -K{y{p) - pb^{vd))- 
M 

Pb is a, known constant, and the dimensionless factor (3 has to be determined experimentally. 
For details see |OTW95a| and |OTW95b] . 

However, a rigorous mathematical derivation of the effective fluid behavior at the boundary 
was not available until Jager and Mikelic applied the theory of homogenization to the problem. 
In |JM96| they developed a mathematical boundary layer together with several corrector terms, 
which allowed them to justify a jump boundary condition. The main tool was the construction 
of several ‘boundary layer functions’: These functions have a given value at the interface and 
decay exponentially outside it. They can be used to correct the influence of spurious terms at the 
boundary, stemming from the contributions of other functions to the fluid velocity and pressure. 

In |JM00| . this theory was applied to give a mathematical proof of the Saffman modification 
of the boundary condition of Beavers and Joseph (see also Section 4 of the Chapter “Homoge¬ 
nization Theory and Applications to Filtration through Porous Media” in |EFM00| for a more 
comprehensible, simplified version of the proofs), yielding the condition 

e (V vpv) ■ T = avp ■ t 


where a = can be calculated explicitely. The constant Cd stems from a boundary layer 

problem for the Stokes equation, cf. [JM00| . Numerical simulations of the boundary layer func¬ 
tions can be found in |JMN01] . 

These results suffer from several drawbacks: First, only a planar boundary in the form of a 
line or a plane is considered (this also applies to the results of Beavers, Joseph and Saffman). 
Therefore, the effect of a possible curvature of the interface is not known. Second, the external 
force on the fluid, appearing as a right hand side in the Navier-Stokes equation, had to be 0. 
This issue was adressed in the recent paper |MCM12] . together with the derivation of the next 
corrector for the pressure. 

Generalizations of the boundary layers in |JM96| were developed by Neuss-Radu in |NR00| . 
However, applications only treat reaction-diffusion systems without flow, and explicit results can 
only be obtained in the case of a layered medium, see |NR01| . 

The main problem which makes the treatment of general settings infeasible is the loss of expo¬ 


nential decay of the boundary layer functions (cf. Section C.2|: With the generalized definition. 


Neuss-Radu was able to show in |NR.00| that an exponential stabilization is not possible in a 
general setting. However, all available tools for the treatment of these problems depend on this 
type of decaj0 

^Maria Neuss-Radu, private communications. 
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Figure 1: Schematic illustration of the velocity profile for a horizontal flow in a domain 
consisting of an impermeable upper boundary (with no-slip condition), a free fluid part 
and a porous region, vp denotes the velocity in the free fluid domain, whereas Ud is the 
effective Darcy velocity. The quantity An = ufIe — vd corresponds to the jump across 
the interface as discussed in Equation 


In this work, we use the approach developed in |Dob09| and |DB10| for providing a gener¬ 
alization of the law of Beavers and Joseph for curved interfaces. We closely follow |MCM12] 
for investigating the effective behaviour of a free fluid in contact with flow in a curved porous 
medium. The main idea is to transform a reference geometry with a straight interface to a 
domain with a curved interface. It is assumed that the porous part in the reference geometry 
consists of a periodic array of scaled reference cells and that the flow in the transformed geom¬ 
etry is governed by the stationary Stokes equation. Therefore, one obtains a set of transformed 
differential equations in the reference configuration. Boundary layer functions for these equa¬ 
tions are constructed such that - due to the straight transformed boundary - their exponential 
decay can be assured. The difference to |NR00| and |NR01| is that in theses works, a periodic 
geometry was intersected by a curved interface. In this work, a periodic geometry with planar 
interface is transformed to give the geometry in which the fluid flow takes place. These results 
have been announced in |Dobl4] and |Dobl5| . In comparison to the latter paper, we will present 
the transformed equations in a more general formulation, hopefully facilitating generalizations. 


2 The Problem on the Microscale 

2.1 Description of the Flow using Coordinate Transformations 

In this section we describe the main geometrical setting which is used throughout this work. Let 
L,K,h > 0. Then := (0,L) x {—K,h) is a rectangular domain in (later corresponding 
to the transformed domain) with parts := (0,L) x (0,h) (later the transformed free fluid 
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Figure 2: Illustration of the coordinate transformation and the periodic domain. 


domain), O 2 := (0, L) x {—K,Q) (the transformed porous medium) and S = (0, L) x {0} (later 
the transformed interface). 

Let g G C°°(U) be a given function such that g(x + L) = g(x) for all x G K. We con¬ 
sider g to describe a periodic curved structure in our domain of interest. Define the coordinate 
transformation 


tjj : fl —> fl 

such that D = V'(^)) ^1 := V'(^i)) ^2 := ' 0 (^ 2 ) and E := '0(S) = {(x,(;(x))|x G (0,L)}. We are 
interested in the behavior of a fluid flowing through the curved channel O, where Di represents 
a domain with a free fluid flow, and O 2 is a porous medium. We are especially interested in the 
behavior of the fluid at the curved boundary E. See Figure]^ for an illustration. 

Let there be given a solid inclusion ^5 CC ^2- We will later use a sequence of such inclusions 
to create a porous medium via homogenization theory. For a given volume force / G L^{Cl), we 
assume that a mathematical description of the fluid is given by the steady state Stokes equation 
with no slip condition on the boundary of the solid inclusion and on the outer walls 


fj,A^u{x) + \7xP{x) = f{x) 

in H\f 25 

( 2 . 1 a) 

diV 2 :(M(x)) = 0 

in H\f 25 

( 2 . 1 b) 

u{x) = 0 

on dCls U 3H\({xi = 0} U {xi = L}) 

( 2 . 1 c) 

u,p are L-periodic in Xi 


( 2 . 1 d) 


Here g > 0 denotes the dynamic viscosity; we will set /r = 1 in the sequel. We are looking for a 
velocity held u G and a pressure p G L^(H)/IR. The Stokes equation is an approximation 

of the full Navier-Stokes equation which is valid for low Reynolds number flows. Using the 
transformation rules for the differential operators (see Appendix 0 . we obtain the following 
equation for the transformed quantities u{z) = u{'ip{z)), p{z) = p{tp{z)) and f{z) = f{tp{z)) in 
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Figure 3: Example of a periodic domain to which our methods can be applied. The upper 
part corresponds to the free fluid domain, and the lower (gray) part to a porous medium. 


the rectangular domain 17: 

-div^(F“^(z)F“^(z) V^u(z)) + F~'^(z) V^ p{z) = f{z) in 17\17s 
diYz{F~^{z)u{z)) = 0 in 17\17s 

u{z) = 0 on 317s U 317\({zi = 0} U {zi = L}) 
u,p are L-periodic in zi 


(2.2a) 

(2.2b) 

(2.2c) 

(2.2d) 


Here 17s •= ‘4’ is the transformed solid inclusion, and F is defined as the Jacobian matrix 

of ■)/) given by 


F{z) 


1 0 
g'izi) 1 


(2.3) 


Since detf’ = 1, ^ is a volume preserving C°°-coordinate transformation. In this connection, 
please note that we define the gradient of a vector field column-wise, i.e. V tt is the transpose of 
the Jacobian matrix of u. Defining V u row-wise leads to slightly different transformed equations. 
Now, the crucial assumption is that 17s is given as an s-periodic structure. This is described in 
the next subsection. 


2.2 The e-periodic Problem 

We assume an e-periodic geometry in 172: Define a reference cell as 

F:= [0,l]^ 

containing a connected open set Ys (corresponding to the solid part of the cell). Its boundary 
dYs is assumed to be of class with dYs n dY = 0. Let 

F* := Y\Ys 

be the fluid part of the reference cell. For given e > 0 such that G N, let x b® the characteristic 
function of Y*, extended by periodicity to the whole K^. Set •= x(^) define the fluid 

part of the porous medium as 


17^ = {a: € 172 I X^{x) = !}■ 
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Figure 4: The reference cell, consisting of the solid part Ys with boundary dYs, and the 
fluid part Y*. 

The fluid domain is then given by 

u s u ni. 

In order to be able to obtain the effective fluid behavior near S, we have to define a number 
of so-called boundary layer problems, see Section [C.2[ To this end, we introduce the following 
setting: We consider the domain [0,1] x M subdivided as follows: 

Z+ = [0,1] X (0,cx)) 

corresponds to the free fluid region, whereas the union of translated reference cells 

z“ = 0{r--("))\s 

is considered to be the void space in the porous part. Here 

5=[0,l]x{0} 

denotes the interface between and Z~. Finally, let 

Z = Z+\JZ- 


and 

Zbl = U S' U Z- 

be the fluid domain without and with interface. For illustrations, the reader is referred to 
Figure]^ Figure]^ and Figure]^ 

Using the constructions given above, we are interested in the limit behavior of the problem: 


For £ > 0, And a velocity and a pressure such that 

-div(F-^U-^ Vu'^) -f = f infl^ (2.4a) 

div(U-iit^) = 0 in (2.4b) 

u® = 0 on 9U®\({a:i = 0} U {a;i = L}) (2.4c) 

u^iP^ are L-periodic in xi (2.4d) 


Note that here and in the sequel, we use the variable x to designate points in U, while y is used 
for the strip Z or the reference cell Y. There exists a solution {u^,p'^) G x L^(U'^)/IR of 

the problem above, see Appendix |C.1| 
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Figure 5: The boundary layer strip ^bl- 


3 The effective behavior of a fluid at a curved porous 
medium 

The effective velocity field in the free fluid domain is given by: Find a velocity u®® and a 
pressure p®® such that 


'^Vu®®) + F-^Vp®® = / 

in Oi 

(3.1a) 

div(F-iu®®) = 0 

in Di 

(3.1b) 

[ p®® dx = 0 


(3.1c) 

u®® = 0 

on (0,L) X {h} 

(3.1d) 

u®®,p®® 

are F-periodic in cci 

(3.1e) 

zi®ff = -eC“ 

on E 

(3.1f) 


C'°'^ is the decay function of the boundary layer function defined in Section 


4.1.2 


It holds 


,eff/ 


The effective Darcy pressure in 0.2 is given by 

div(F-iy4(/ - Vp®®)) = 0 
A{f- F-'^ V p®®) • F-^e2 = 0 


~eff 


direction, a slip velocity 
n Lemma |C.23| as 

-£(7“ • Fei 

d^\x,yi,+0) ■ F{x)ei dpi 

(3.2) 

in 0,2 

(3.3a) 

on (0,F) X {-K} 

(3.3b) 

on E 

(3.3c) 

is L-periodic in xi 

(3.3d) 









A is given in (4.5), and C^* is the pressure stabilization function defined in Section 
by 

= [ w“(a:i, 2 /i,+ 0 ) dyi. 

Jo 

Define the effective mass flow rates in transformed tangential direction as 


4.1.2 


given 


:= / • Feu dx := / u" • Fei dx. 

We obtain the following estimates: 

3.1 Theorem. 

Let f € be L-periodic in the first variable. For as defined in (2.4) and 

defined in ([33 the estimates 


2 + \M‘^ - < Ce^ 

I 11 T (f2i) 

^efl 


M® —u®®|| 1 o + \\p^—p^^\\ 1 + ~ w®“)|| 1 i<Ce 

rJ-o,•n.^2 V / r.d/'Onl^ — 


eff N 


'H 2 (f2i) 




|a;2|5 




\x2\Hf-Pn 


< Ce 




hold. 


In the porous medium 122, we arrive at the following results: 

3.2 Theorem. 


For the effective pressure in the porous medium defined by (3.3), we have for all S > 0 
- A{f - F-^ Vp^^) —^ 0 m L^{{0,L) x {-H,-S)) 


3.3 Theorem. 

At the interface E, it holds 


P -P 


~eff 


p -p 


•eff I 


0 in L^{VL 2 ) 
, _i <Ce5 

Ijf 2 (E) — 


(3.4) 




,.eff I 


Ih 2(2) 


< Ce-^. 


(3.5) 


These results show that the following behavior of a free fluid in contact with a flow in a 
curved porous medium can be expected for low Reynolds number flows: 


• In the free fluid domain, the velocity and pressure are given by the Stokes equation. 

• In the porous medium, the flow is pressure driven and given by Darcy’s law. 

• At the interface, a slip-condition occurs. The velocity is given with the help of the decay 
function C^^{x) of an auxiliary boundary layer problem. This depends on parameters from 
the geometry of the interface. 
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In transformed normal direction C^^{x) ■ F~^{x)e 2 = 0 holds, which is an approximation 
of continnity of the velocity in that direction. In tangential direction, a .jump between 
the velocities and pressnres occnrs. Looking at the form of the bonndary condition, we 
see that the generalized bonndary condition of Beavers and Joseph has to incorporate 
effects stemming from the geometry of the flnid-porons interface. The estimates (3.4) 
and (3.5) corresponds to a generalized rigorons version of (1.2). In |JM00| and |MCM12] . 

the jnmp bonndary condition is given by uf® + = 0, whereas in onr case it 

reads u®® + = 0. The missing velocity gradient appears in the eqnations for see 

Section 14.1.21 


4 Derivation of the general law of Beavers and Joseph 


In this section we carry ont the steps that are necessary to derive the generalized bonndary 


condition of Beavers and Joseph (3.2) as well as the theorems given above. We will snccessively 
correct the velocity and the pressure p® given by (2.4) with the help of anxiliary fnnctions. 


This will give ns insight into the effective behavior, while at the same time introducing problems 
with div(i^“^M®). Therefore we have to correct this term as well. 


4.1 Correction of the velocity 

4.1.1 Elimination of the forces 


We start by eliminating the right hand side of (2.4a) in fli. Let be a solution of 




Vu°)+F-^Vtt° = / in fli 
diY{F-^u°) = 0 in fli 

= 0 on i9ni\({a;i = 0} U {xi 
vP,Tr^ are L-periodic in xi 


L}) 


(4.1a) 
(4.1b) 
(4.1c) 
(4. Id) 


There exists a unique solution G 7r° G L^(ni/K) by the results for the transfomed 

Stokes equation. By regularity results (see e.g. |Tem77| L this solution is smooth for smooth /. 
We extend the velocity by 0 in ^2 and the pressure in a smooth manner to a pressure 7r° 
defined in all of 11. As it will turn out, tt® will be given by the Darcy pressure in D 2 . Details of 
this extension procedure will be given below (see Section 4.4). 

To obtain estimates for — u^, we need the following Lemma from |JM96| . which is a variant 

of the Poincare inequality: 


4.1 Lemma. 

Let (j) G {g G FI^{n 2 )\g = 0 on 9D2\i9D2}. Then 

as well as < Ce ||V ■ 

Define the space = {g G i7^(n®)^ \ g = 0 on 90®\({xi = 0} U {xi = 

L}),p is L-periodic in xi}; then by subtraction the weak formulations, one easily sees that 
— vP satisfies the variational equation 


/ F ^(V — V M^) : A ^ V (j) dx — f (p® — 7r°) div(F dx 


= / F-^F-'^Vu° 

JY. 


e2-4>da^- f F • 62 dcr^, + f {f - F V Tt°)4> dx 

Jy Jog 


(4.2) 
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for all (j) G W^. Here [g]s := —q\a 2 denotes the jump of the function q across the boundary S. 

4.2 Proposition. 

We have the estimate 


y/e\\Vu^ -Vu° 








Proof. Choose cj) = in (4.2). Since for u,v G H^{VP) the form a{u,v) := F ^ Vm : 

F~'^ V n da; is bounded and coercive, we obtain 


Vu® - Vit' 


0| 


I 


<C f |F-^(Vm'^-Vm°)P 


da; 


< 

[ F ^F ^ \7 it°e2 • M ~ do’s 

+ 

f F ^[7r°]E(lt® - lt°) • 62 dtTa; 


JY. 


Jy 


[ da 

Jm 


Since V and are smooth, the first two terms on the right hand side are bounded by 


C'||m'^-u°|| „ , < Cv^lIVu" - Vu° 




The function {f — F V 7r°) is smooth as well by elliptic regularity results for the Darcy pressure, 
so the last term is bounded by 


CM -u^ 


K , 2 < Ce 

\L^{n^y — M 


I L^{Q^y 


< C{S)e^ + CS\\Vu'^ -\/u° 




Choosing S small enough we arrive at 11V — V 
1, 


= C£ II -Vu° 


< Ce. Moreover, we have 
< Ce^ 


2/'o?^4 — 


I 


as well as 


IIm''||^ 2,^.2 < Ce ||Vu® - Vito||?2,„e>4 < Ce^. 






This finishes the proof. 


□ 


In order to derive estimates for the pressure and it® — in Hi, we use the theory of very 
weak solutions for the transformed Stokes equations, see Appendix [B| 


4.3 Lemma. 

It holds 


/illp® —7r°|| , +||it® —it^ll , , - 


< Ce. 


Proof. It® — is a very weak solution with Gi = G 2 = 0, .^ = '«®|e,j,, St = L U {a 2 = h}. Using 
the estimate (B.4| from the theory of very weak solutions of the transformed Stokes equation, 
we obtain 

<C'||eiL2(,^)=G|k®|L2,.. <Ge 


\u^ — u 




IL^(S) — 


by the previous Proposition. Using the Necas inequality yields 
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4.1.2 Continuity of the traces 


Looking at the right hand side of equation (4.21 and the proof of Proposition [4^ one can see that 
the expression — V 7r°)(/) da; allows estimates for ||Vit® ~ ^on the order ofe, 


whereas the other two integrals only allow for estimates on the order of eJ. Later we will choose 
to allow for better estimates. This leaves us with the expression V u°e2 • 4> dax, 

which we are going to correct next: Construct the following parameter dependent boundary layer 
functions satisfying 


-divy(F ^{x)F '^{x)'Vy P'°\x,y)) + F (x) Wy ui'^\x, y) = 0 

in 0 X Z 

divy{F-'^{x)/3'^\x,y)) = 0 

in H X Z 

[/3“(x,2/)]s = 0 

on H X S' 

[{F-\x)F-^{x)\'y/3'^\x,y) - F-^ {x)uj^\x,y))e2]s 


= F~^{x)F~'^{x) V u°{x)e 2 

on H X S 

0^\x,y) = Q 

on H X Ur=i{5>^S - (Dl 

0°^{x, •),w'’*(a;, •) are 1-periodic in yi 



and define = e/3*’^(x, |) as well as = a;’^'(a:, |). 

By the theory of the boundary layer functions (see Appendix |C.2[ ), there exists decay functions 

Cbl . ]^2 -^ : m2 -^ 


I ||/3“’^-eC'“(x)iL(x2)|L,„^2 + -C'“(a;)iL(a;2) 


V,(/3“(a;, ^) - iL(x2)C'“(a;)) 
£ 


4- 




I L^(n) 
3bl/ 


|Vy/3“’"| 


div,(F-^l^-^ VyP^^x,-)) 


< Cel 




Here H denotes the Heaviside function. This correction introduces problems due to the stabi¬ 
lization towards Therefore we define the following counterfiow: 


— div(A ^{x)F ^(x) V m'^(x))- 1-S' ^(x) V 7r°’(x) = 0 

in Oi 

div(F“^(x)u°’(x)) = 0 

in Hi 

e 

<5 

II 

0 

on (0, L) X {h} 

U'^{x) = c'^\x) 

on E 

u'^,Tr'^ are L-periodic in Xi 



Since J^C'^^(x) ■ F~'^(x)e 2 dx = 0 (see Lemma C.21I, there exists a unique velocity u'^ and a 
pressure unique up to constants. Consider H{X 2 )) + eu'^H{x 2 ) as a first 

macroscopic approximation. Define the error between this approximation and the microscopic 
velocity (similarly for the pressure) as 


W = u^-u° + (/3“>^ - eC'“iL(a;2)) -f eu'^H{x2), 

=p^- tt^H{x 2) - n°H{-X2) + (w“’" - C'“iL(;c2)) + eTT'^H{x2). 
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Then we have 


/ F \7U^ : F ^ (j) dx — f div{F ^(f)) dx 

Vne 

= - f T'“^[^°]Ee 2 • (/> dcr - / F~^C^^e 2 ■ 4> da + f {f - F~’^ V dx 

Jn Jy. Jai 

+ [ -eC'^'H{x 2 )):F-^V(l)dx- f div^{F-^F-^ Vy p'^'ix, -))(!> dx 

Jqs £ 

+ e [ F-^F-'^ V u'^e 2 ■ p da + [ - F(x 2 )C'“) • </> dec + e [ F-\'^e 2 ■ p da 

Jy dn® Jy 

+ [ F-iF-^V,(/3“’"-eC'“F(ec2))e2-<^d(T 

J {X2=h} 

- [ F-iF-^V,(/3“'"-eC“F(ec2))e2 da. 

J {x2=-H} 


We will later eliminate the first two terms on the right hand side by requiring that = —C}^ 
(denoted as (P), see below). Due to the exponential decay of the boundary layer functions, the 
last two terms can be chosen arbitrarily small (denoted (BL)). The terms on the right hand side 
can be estimated by the following orders of e with respect to VU’^ (the signs of the terms are 
kept in order to facilitate the allocation of the terms, and \\ denotes a line break): 


F-'^\7U 


+ 


: F ^ \7 (j) dx — j div{F ^0) da; = — (P) — (P) + i^(e:) 
dO' 

ei)-^(£5) \\ +^(£i) + ^(ei) + ^(£i) \\ +(bl) \\ - (BL). 


Now is of order es, and is of order £3 in For the proof of the main result, we need 
to be on the order of £, and W should be on the order of . We therefore need to correct 
the velocity again. 


4.1.3 Second correction of the velocity 

Define the following boundary layer function: 


-diVy(F ^(a;)F {x)Vy')^\x,y)) + F {x)^y y}^\x,y) 

= diY,,{F~^[x)F~'^{x)\Jy /3“(x,y)) 


div,(F-i(x)7“(a;,2/))=0 

[7“(2;,y)]s = 0 

[{F-^{x){F-'^{x) Vy 7“(x, y) - i?\x, y)))e2]s = 0 

l°\x,y) = 0 

7'^^(a:, •),/i’^*(a;, •) are 1-periodic in yi 


vD. Q, X Z 
on fl X S' 
on D X S 
on n X UJli 


Ors-ffl} 


in Q X Z 
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The decay functions are denoted by Define the corresponding counterflow 


— div(F ^(x)F" ^(x)V6(x)) + F '^{x)Vq{x) = 0 

in Di 

div(F“^(x)6(x)) = 0 

in Oi 

II 

o 

on (0, L) X {/i} 

6(x)=C'“(x) 

on E 

b, q are L-periodic in Xi 



4.1.4 Correction of the pressure 

For the correction of the pressure, define the following boundary layer fnnction 


-divy(F ^{x)F '^(x) Vj, A“(x,?/)) + F ^(x)VyK“ 

{x,y) 

= F-^(x) V,(c.“(x,y) - iJ(x 2 )C'“(x)) 

in D X Z 

diVy(F"^(x)A“(x,j/)) = 0 

in D X Z 

[A“(x,2/)]s = 0 

on D X 5" 

[(F"i(x)(F"^(x) Wy A“(x, y) - k“(x, y)))e 2 ]s = 0 

on D X S' 

A“(x,y)=0 

on D X Ur=i{'9>S- (°)} 

A’^*(x, •), k’^'(x, •) are 1-periodic in yi 



The decay functions are denoted by Define the corresponding connterflow 


— div(F ^(x)F ^(x) V c?(x)) + F ^{x)VZ(x)=0 

in f2i 

div(F“^(x)d(x)) = 0 

in 111 

c?(x) = 0 

on (0, L) X {h} 

d{x) = C\\x) 

on E 

d, 1 are F-periodic in xi 



We extend onr approximation to and by adjnsting lA^ and . Define 

W=W + - e"i?(x2)C“) - *) - e^H{x2)C'f) + e^H{x2)b - e^H{x2)d 

8 ' 8 

= + {£p'°\x, - eiJ(x2)C“) - (eAc“(a;, - eH{x2)C^^) + e^H{x2)q - e^H{x2)l 

The weak formulation is now given by 

/ F~^ W ■. F~'^ (j) dx — f div{F~^(j)) dx 

= - f F~^[TT°]^e 2 ■ (j) dcr - f F“^C“e 2 • <^ dcr + f (/- V 7f°)(/) dx 
Jn Js Jn^ 

+ [ F-^ - eC“71(x2)) : V dx + e [ F-^T:'^e 2 ■ </< do- 

+ £ I F~^F~^ \7 u'^e 2 ■ (jy da 

Jt, 
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+ [ e2i^-^V,(7“(x,-)-i7(:r2)C“) :F-^V(^dx- [ e div,(F-iF-^ V,, 7 “(x,-))(^ dx 
dn® £ Jn= ^ 

+ [ e‘^F-^\7^{X'^\x,-)-H{x2)C^'):F-^V<l)dx-f ediv^(F-iVy A“(x,-))(^ dx 
e Jns e 

- / F-^F-'^ Vib - d)e 2 ■ (j) da + / VJfJ^Ux, -) - H(x 2 )C]!^U dx 

ds Jne e 

+ [ eF~^C^e 2 ■ (j) da — [ eF~'^ — H{x 2 )C]’^)(j) dx 

dE dne e 

— f eF~^C]^^e 2 ■ (j) da + f — Z)e 2 • 0 dx + (B). 

JE JE 

Here (B) designates terms stemming from the boundary layer functions at the outer boundary 
of H. Again, due to their exponential decay, they can be chosen arbitrarily small. Similar to the 
calculations above, we arrive at the following estimates of V in terms of e: 


F-^WW : F-^ Web dx- 


-T . 


V‘^ div(F ^(j)) dx = —(P) — (P) + ^(e) 


+ ^(£S) + ^(eS) \\ +^(eS) \\ +^(e5)-^(£S) \\ +^(£S)-^(£S) 

_^(£i) + ^(£f) \\ +^(ei)_^(ef) \\ -^(£i) + ^(£S) + (BL). 


Thus, one arrives at an approximation ~ in HI and W\y. ~ In 

order to prove these estimates rigorously, we need to correct the divergence of U’^ next. 


4.2 Correction of the divergence 

4.2.1 Compressibility due to the boundary layer functions 

A calculation shows that 

div(A-iW^)(x) = £div,[A-i(/3“(x, - it(x 2 )C'“(x))] + £2 div,[F-i( 7 “(x, - H{x 2 )C^\x))\ 

e e 

-e^ diva,[J^“^(A“(x, -) - H{x 2 )C^\x))] 

£ 

We construct the following convergence of the divergence: Let satisfy 


diVj^(F ^(x)(3/j(x,?/)) = diVa,(^J^ ^(x) [/3“(x, y) - iL( 2 / 2 )C“(x)]) 

iwVtxZ 

Q/3{x,y) = 0 

onHxUr=i{^5-(D} 

[Q 0 ]s{x,y) = Cf{x) 

on H X S' 

Qp{x,y) is 1-periodic in yi 



where 

Csix) = F{x) [ div,j,(A"^(x)(/3“(x,?/) - H{y 2 )C^\x)))e 2 dy. 

•^Zbl ^ ^ 
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Set Qp{x) := e^Qpix, |). Similarly, we define <37 and Q\ 


divy(F ^{x)Qj{x,y)) = diY:,(^F ^{x)['y'^\x,y) - H{y2)C'^\x)]^ 

in n X Z 

Q'y{x,y) = 0 

onnx[jZi{ys-0} 

[Q.y]six,y) = C^{x) 

on fl X S 

Q^{x,y) is 1-periodic in yi 



as well as 


divy{F ^{x)Qx{x,y)) = div:,(^F ^{x)[X'^\x,y) - H{y 2 )C'^\x)]'^ 

in fl X Z 

Q\{x,y) = 0 

onOxUr=i{^S-(D} 

[QA]s(a:,y) = C^{x) 

on 0 X 5” 

Qxix,y) is 1-periodic in yi 



with decay functions 


C^{x) = F{x) / ^{x){-f'^\x,y) - H{y 2 )C^\x))'^e 2 dy 

Cxix) = F{x) f diY^(^F-^{x){X'^\x,y) - H{y 2 )Cl\x))^e 2 dy. 

Again, define Q^(x) = e^Q^(x, |) and Qx[x) = e^Q\{x, |). The tools and techniques for proving 
existence and exponential decay for these type of problems can be found in Appendix |C.3| 

The jump across the boundary of the functions above is corrected with the help of the following 
counterfiows: For i € {/ 3 , 7 , A} define as the solution of 


— div(F" ^{x)F '^{x)'\/u^'{x)) + F ^(cc) V (x) = 0 

in rii 

div(F-i(a;)u'3^(^)) = 0 

in Oi 

v9'{x) = 0 

on (0, L) X {h} 

u^'-{x) = Cf[x) 

on S 

yQi^jj-Qi are L-periodic in xi 



These problems have a solution if Cf ■ F ^62 dcr = 0. For we obtain 


C9 ■ F-^e 2 da = 


f f div:j(j^ ^{x){p'^\x,y) - H{y2)C'^\x))) dy dao; 
j diVa:(F"^(a ;)[J f3'^\x,y) - H{y2)C'^\x) dy]^ da^; 

/3“(^,y)-ff(2/2)C'“(x) dy]) -ei 

+ J - H{y 2 )C'°\x) dyj) • 62 da^; 


Now the first term on the right hand side vanishes due to the periodic boundary funcion y, 
which leads to periodic boundary conditions in x for F, and The second term vanishes. 
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since the latter functions do not depend on X 2 - The argument for the remaining counterflows is 
analogous. 

Now define 

Ul = + £^vP-< - £^v9^ 

7?^ = -p® + £2^Q/3 

which leads to 

div(J^“^(x)W®(x)) = -e^ AW:^{F~^{x)Qp{x, -))-£^ diva;{F~^{x)Q^{x, -))+e^ div:,{F~^{x)Qx{x, -)) 

£ £ £ 


4.2.2 Compressibility due to the auxiliary functions 

We correct the divergence even further by defining the functions and via 


div(F ^(x)(()-’®(a;)) =-Ei diva;(F ^{x)Qi{x,^)) 

in Hi 

o 

II 

on (0, L) X {h} 

{x) - [ £idiv^{F ^{x)Qi{x, ))e 2 dx 

1^1 Jni ^ 

on S 

is L-periodic in xi 



and 


div(i^ ^(x)0-’^(x)) =-E 4 divj;(P ^(x)Qi(a:, ^)) 

in H 2 

P1’’^{x) = 0 

on (0,L) X {-K} 

4>i’ [x] - I I [ £idivyF ^{x)Qfx, ))e 2 dx 

1^1 J 02 ^ 

on S 

is L-periodic in xi 



Here Qi is extended by 0 into ~ (D}> * ^ and 

{ e^ for i = /3 
£^ for i = 7 . 

—for i = A 

It holds < C\£i\, < c\£,\ and 4>Y = -</>■’" on S. 

We define the final correction by making use of the following restriction operator (see |Dob09| 
for details): 

4.4 Proposition. 

Let H = {u & : u = 0 on {0, L) X {—F[},u is L periodic in the first variable} and set 

Fl^ = {u & H ■. u = d on i9H2\5f22}. There exists a linear restriction operator TiT : FI —> 
such that for w G H: 

w = 0 on dLl^\dn TZ'^w = wine 

div(F“^w) = 0 in fl =7 div(P“^7?,^'u;) =d inLF 


and 


iTZ'^wl 




■e||V(7^^w) 




4 < Cdlwl 


L^(ny 


+ £ ||v w| 
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Using an explicit characterization of 7?,® (see the Appendix of |SP80| by Luc Tartar), one 
arrives at the following identity 

= div(F"V) + f div(F"V) di. (4.3) 

e(>s+fe)cn2} 


Now define 


W^=UI-H{x2) Y. - H{-X2) Y 


2 ,£ 


n = n- 


Since diva;(F ^[x)Qi{x, |)) = 0 in e{Ys + fc), /c G 1?, equation (4.31 leads to div(U ^lA^) = 0, 
and we can use Uq as a test function. We look at the weak formulation 

f F~^ \7 Uq ■. F~^ \/(f) dx — ( T^q div(U“^0) dx. 

JOS 

By calculating the right hand side of this equality and inserting (j) = Uq , we obtain (similar to 
the proof of Proposition 4.2| < Ce, < Ce^ and \\UI,\\^ 2 < Cei. 

Similarly, one arrives at ||'Ug (lif-qaC) < Ce. 

4.5 Corollary. 

Let 7r° he a smooth function satisfying = —C^^. Then the estimates 


\\W\\ 1 ,+e--\\V^ 

" ff 2 " 




hold. 


Proof. The result follows as in |MCM12] by transforming the equations in Ui back to a Stokes 
system in Ui. □ 


4.3 Some Estimates for the Corrected Velocity and Pressure 

Note that differs from Uq by ^(e^) in the L^-norm, and by ^{e) in the iL^-norm. This yields 
the estimates 


ellVP" 


qs7^)' 


+ e||VW" 




\\W 


\L^(a^y 


+ £= \\w 




as well as 


e\\r^ 




\\u^ 


' H 2 (Qj^) 


,<CeF 


We now collect some more estimates. They are analogous to what is known in the literature 
|MCM12] . and the proofs follow along the same lines. 


4.6 Lemma. 


\u — u 


eff I 




< Ces 
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Proof. Note that — u®® = W — — eH{x 2 )C'^^). By the theory of very weak solutions we 

have 


I lli,2(ni)2 


< C \\W - (/?“’" - < CeS. 


□ 


4.7 Lemma. 

Proof. First note that 




< l|VZ^1L2,,^,2 + ||V(/3“--£i7(x2)C 


'bl\ 


Ii^(r2i)^ _ I, • -- Iii,2(ni)^ 

The second assertion follows via 


2 fo, \2 — 




< Ce + Ce3 < Cek 


||V(n®- 

4.8 Lemma. 


s „,eff 


U 


iL^iniY 


< I1VW®|L2(,^,2 + ||V(/3“’® - £i7(x2)C'“)k,^ , < Ce. □ 


|u® — u®®|| 1 , < Ce 

I lljf 2 (f2i)^ — 


Proof. By interpolation we obtain 


u® - M®®|| 1 <C Ik® - U 


€ ^,eff II 2 


Ih2 (S2i)2 




V(m® - M®“)|| ^ 2,0 ,2 < Ce^'^ei'i < Ce. □ 




4.9 Lemma. 


k2kVk®-U®“) + k2kV(p®-p®«) 




< Ce 


i^(S2l) 


Proof. We have 
k2|5 Vk®-U®®) 


< Ce + 


< Ce -\- C ||Z^Q II + |X 2 | 2 


ksl^ VU^ + |x2|i - eH{x2)C^') 






L°°(Qi) 




For the second part, nse Proposition 16 in |MCM12] to obtain 


\X2\^ V(p® -p®®) 


< c|k^-w®®' 






+ Ce < Ce 


□ 


4.10 Lemma. 


Proof. 


\\p" 


Ik' ^ 11^' - - HiX2)C^') - eiFk2)^kLi 




< imiL^tn,) + - ^^k2)C^k.p.p + Ce < Ce 


□ 


4.11 Lemma. 


|M® - M®®| < Ce* 
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Proof. 


4.12 Lemma. 


\M‘^ 


= {u^ - ■ Fei dx\ 

< C IIm® — u®®|| 2 ,< Ce^ 

— II Ili2(s2^) — 


□ 


1 

e 





\u — u 


eff I 


—^ 0 m L^{Yf, 

I _i <C'e5. 

' H 2 (S) — 


Proof. It holds — u®® = W — — eC'^*) on E. Since ||£ 2 , 2 ( 5 ,) ^ Ce and -) —^ 

Jo 2 /I 1 +0) dt/i = C'’^*(-) in L^(S) (cf. |Amo97| or the remarks in |A1192 | ]. we obtain the first 
result. 

8 

Now let V S H^(T,). To begin with, we have 




= e 


Qh\ 


iP (syij+0jyii+0) 


(/ 3 “(xi, +0, —, +0) - c“(xi, + 0 ))?;(a;i) dxi 

£ 

- C'“(eyi,+0))z;(eyi) dyi 


Jo 

= ^Y] [ {P'°'{eyi,+0,yi,+0) - C'^\£yi,+0))v{£yi) dyi, 


where 7® C No is some finite index set depending on e, and Yi is an interval of the form \k, fc + 1) 
for some k e Nq. For each i G 7®, choose ay* G Since /3'^'(y*, +0, •, +0) — +0) has mean 

value 0 over E, we have that {[}'°^{£y*, +0,yi, +0) — C'’^'(ey*, +0)) • v^ey *) dyi = 0. Thus the 
expression above is equal to 

eV / (^“(eyi,+0,yi,+0) - C“(eyi,+0))r;(eyi) 

^£Ie Jyi 

-(/3“(£2/z*>+0,2/i,+0) - C'^\ey*,+Q))v{£y*) dyi. 


By Taylor series expansion, there exists a ^* G Yf such that the above expression is equal to 


^ X! / +0’ yi> +0) - +0))^^(eOj e(2/i - vt) dyi 

,e/e JY! oxi L J 

= eX! 4,^ (/3 “(£^Cz*)+0:2/i>+ 0)-C'“(eC,+0))ii(e^r) m dyi 


iG/' 


by the vanishing mean value of /S'^* — C^*. Now transformation of the integral and rearrangement 
yields 


I f7^[/^“(4*,+0,-,+0)-C'“(e*,+0)]i;(er) 

t^eJeYl'-OXi e 


+0, y, +0) - +0) 


Xi dxi 
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By the embedding ^ C^(K) as well as by the boundedness of the terms containing 

^bi _ integral can be estimated by eC ^ which upon summation gives 

a bound of eC||v|| s for the whole term. This implies that — eC'^'ll _8 < Ce^. 

" "ff5{i:) ^ Ir llff 5 (s) — 

Since ^ we have that ^ Ce^. Thus similar to the first part 

||y£ _ yeffll Ce^, and by interpolation we obtain 

II II rr 5 (s) — 


.eff I 


Iff 2(S) 


<C 


,eff I 




.eff I 


3 

■2 (S) 


< = C£5. 


□ 


4.4 Proof of Theorem 13.21 

We need the following auxiliary constructions: Fix a; € fl and let {w^jTt^) G x L^(V*)/R 

be a solution of the parameter dependent cell problem 

(4.4a) 
(4.4b) 
(4.4c) 
(4.4d) 

Define the matrix A by 


— divy(F“^(a;)F“^(x) Vy w^(x, y)) + F~’^{x) Vy Tr^{x, y) = ej in Y* 

d.iYy{F~^{x)w^{x,y)) = 0 in Y* 

w^{x, y) = 0 in Ys 

{x,y),Tr^ {x,y) are T-periodic in y. 


[A{x)\ji = wl{x,y) dy i,j = 1,2. (4.5) 

One can show that A is symmetric and positive definite. The Darcy pressure is given by the 
following problem: Find 7r° G F[^{VL 2 ) such that 


div(F-M(/ - F-'^ V it°)) = 0 

in Q 2 

(4.6a) 

4(/-F-^Vit°)-F-^e2 =0 

on (0,L) X {-K} 

(4.6b) 

0 

II 

0 

+ 

on E 

(4.6c) 


is L-periodic in xi 

(4.6d) 


Here is given in (4.11. Using the weak formulation for W,V^, one arrives at the estimates 


\\W 




< OeU 


|VW^|U2(^e,4 <C£, 


\\V^ 


|l,2(sj|) < C. 


By the theory of two scale convergence (see |Ngu89|, IA1192n there exist limits iY*™ G 
L^(fl 2 ;B#(y*)) as well as U'™ e ^^(y*)) with 

S 

Choose a test function cj) G C((°(H 2 ; F[^{Yf)Y with divj,(F“^(a:)(()) = 0 and test the weak formu¬ 
lation for . In the 2-scale limit, one obtains 


F ^(s) VyZ^‘™(a:, y) : F \x)yy (l){x,y) dy dx 


ylim 


-T/ 


/fio ^Yf 


[ [ T’'™(x,y)div,i,(F ^{x)(t){x,y)) dy dx = [ I 

JEI2 ^Yfr JElo J 
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I (fix) - F '^{x)V TT°{x))(l}{x,y) dy dx. 

Yf 

(4.7) 











By choosing a test function (p € {0,2] H^{Yf)Y, testing with £(p{x, |) and taking the 2-scale 

limit, one arrives at V^^'^{x,y) divy{F ^{x)(j){x,y)) = 0, thus 7^'™ does not depend on 

the variable y. Upon a separation of scales similar to the derivation of Darcy’s law (see |Hor97| 
and |Dob09| for the case of the transformed equations), one obtains the representation 


U 


““(x, y)=Y. w^{x, y) [fix) - F-^{x)w(n°{x) + 


i=i 


U"’^ix,y) dy = A{x) /(x) - F"" (x) V -k 


Choosing as a test function yields < Cei and thus V ^ —>■ 0 in 

Choose a. (p € {O; H{^{Yf))'^ such that divy{F~^{x)p) = 0. Similar to the above calculations, 

testing with p{x, |) gives in the limit 


/Qa JYf 


F ^ {x)VyW"^{x,y) : F ^ {x) Vy (p{x, y) dy dx 


- f f U'™(x)div,j,(F ^{x)p{x,y)) dy dx = [ [ {f{x)-F ^(x) V 7r°(x))(/)(x, y) dy dx. 

Jn JYf Jq.2 JYf 

Using (4.7l, one arrives at Jy^ V^'™{x) divx{F~^{x)p{x,y)) dy dx = 0. Upon an integration 
by parts, this is equivalent to 


/Sli JYf 


F {x)W''™{x)p{x,y)dydx+ / / V''™{x)p{x,y) ■ F (x)e 2 = 0. 


E JYf 


Since (p is arbitrary, we conclude that VU'™ = 0 in Oi and F'™|e = 0. Thus F*™ = 0 and 
Jy^ Z//'™(x, y) dy = A{x){f{x) — F“^(x) V 7r°(x)), which is Darcy’s law. 


It remains to compare with p®”. Note that 






eff ^01 








After transformation to Ui, the functions it®”, p®”, and p^ satisfy a Stokes equation. By stan¬ 
dard regularity results for this equation (see |Tem77| . |MCM12] 1 we obtain lip®® — tt^II ,, < 

Ce for all k gN. This gives the first statement of the theorem. Additionally, since 

V^ln, =p®-ir° + ^(e5) 

= p® - p®® -f i^(e5) —" 0, 

we have that p® —^ p®® in F^(U 2 ), which is the second assertion. 

Finally, we estimate ||p® — 7r°|| Let w € 77^(171) such that i/; = 0 on 9U2 \({xi = 

0}U{xi=F}UE) and w is F-periodic in xi. Then by trace estimates |Tem77| 

For the first term on the right hand side, notice that p'^ — = (p® — p®®) + (p®® — 7r°) = 

^(eJ) -I- ^(e). For the second term, observe that 

^ ... +|(V(.;“’®-F(x2 )C“),zr;),2(,e,| 


|(p®-7r^w|E)^_l ^l| < |((p®-it^)e2,Vw)^2(s,e)| + |( —(p®-it'J),w)^2((,e)| 


[-^{p" - ^ ),w)F 2 ^^eJ < C'e||rt;||„i,„e, + || VF®|| IIWII 


I(n|) 


I II ^ II 


< C'e lkLi{s2i) +Ce 5 . e IIV w||^2(„e, ■ 


This shows that p® — 7r° < C^/e, which implies the last statement of the theorem. 
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A Coordinate Transformations of Differential Operators 

We recall the definition of coordinate transformations and some differential operators and inves¬ 
tigate their relations: Let il C M" with n G N be a Lipschitz domain; let c : n —)- M be a scalar 
function, j : (l —K" a vector field and M : Ct —a matrix function. They are assumed 
to be sufficiently smooth. 

A.l Definition. 

The gradient of a veetor field is defined as 


^ djk 
dxi' 

for i,k = l,...,n (i.e. V j is the transpose of the Jaeohian matrix off); the divergence of a 
matrix-valued function is defined column-wise, thus 


(Vj) 


(div(M))fc = 

2=1 


d^k 

dx^ 


for k = 1,... ,n; and the Laplacian of a vector field is given by 

Aj = div(Vj). 


For n = 2 we define the two operators 


Curl(g) 



Vc 


and 

dj2 _ ^ 

dxi dx2 


curl(j) 


The ‘curl’-operators above are two-dimensional variants of the well-known curl operator de¬ 
scribing the rotation of three-dimensional vector fields. We have the relations 


curl V c = 0 and div Curl c = 0, 

and curl is the formal adjoint of Curl (see |Vern7| and |DL90| for details concerning these oper¬ 
ators). 

A.2 Definition. 

Let n, C M” be Lipschitz domains and let tp : fl —> Ll. We call ip a regular orientation¬ 
preserving C^-coordinate transformation if 

1. Ip is a -diffeomorphism, and 

2. There exist c, C > 0 such that 


c < det F{z) < C Wz £ LI, 
where F denotes the Jacobian matrix of ip. 

If det F = 1, we call ip volume preserving. 
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We will indicate coordinates in by z = (zi ,..., z„) and those in by a: = (xi,..., a:„). 
Define 


c{z) := cii^iz)) 
j{z) :=j{ip{z)) 
M{z) :=M(V'(z)) 


A.3 Lemma. 

Let 'll) : fl —> LI be a -coordmate tra-nsformation. De'note by F the Jacobian matrix of 'll), and 
let J{z) := det(-F'(z)). With the notations and definitions above it holds 

1. F-'^y^c= VajC. 

2. AFr^{JF-^j) = (Jo'0-i)diva;(j). 

3. z{J F~^ M) = {J o'il)~^) dW x{M). 

Proof. The first assertion is a simple application of the chain rule, whereas the second one is 
known as the Piola-transformation (see |Zei88| . Chapter 61. Note that Zeidler defines vectors 
and gradients row-wise, leading to slightly different formulas.) For the matrix divergence the 
second statement holds column-wise. □ 


Application of this lemma yields: 


A.4 Lemma. 

Let 'ip be a volume-preserving -coordinate transformation, 
transform according to 


The operators from Definition 


A.l 


1. 

^x(c) = 

divz{F-^F-^Wzc) 

2. 

II 

<] 

diYz{F-^F-^\/zj) 

3. 

diVa;(j) = 

= divz{F~^j)- 

1 

divx(M) 

= divz{F-^M). 

5. 

Curl2.(c) 

= Curl;,(c), 


with 



Curl2(c) 


0 -1 

1 0 




6. curU(j) = curl 2 (j), 

'with 

curU(j) = curU(i^^j). 


Proof. For volume-preserving coordinate transformations it holds J = 1, thus in that case by 
the preceding lemma we have div 2 (F“^j) = diva;(j) and d\Vz{F~^M) = diVa;(M), which gives 
the third and the fourth statement. The first and the second statement follow by the equalities 
Ai-xic) = diva,(Va;c) and Ax{j) = div 2 ;(V 3 ;j) and application of the above results to the right 
hand sides. 

The fifth statement follows along the same lines, whereas the sixth can be obtained by a 
direct calculation of the effect of the transformation on the defining equation. □ 
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A simple computation shows that 


; thus it holds 


0 

1 


-1 

0 


F-t 



-1 

0 


Curl^(c) = FCurl 2 (c). 

A.5 Lemma (Transformed Differential Identities). 

Let Ip be a volume-preserving -coordinate transformation as above. Then the following identities 
hold: 


1. div2(F ^c) = F '^VzC. 

2 . div,(T’-ii^-^V,(div,(F-ij))) =div,(F-idiv,(F-iF-^V,j)). 

3. diYz{F~^ Curl^(c)) = 0. 

4- diYz{F~^{cj)) = cdivz{F~^j) + F~'^ V c- j. 

5. curl 2 (A“^ Vz c) = 0. 

6 . curl2(div^(F“iF“^ j)) = divz(A“iF“^ curl2(j)). 

7. If diYz{F~^j) = d, then 


F Vz{cur\z{j)) 


0 

1 


div,(F-iF-^ \/zj). 


8. F-’^Wz{diYziF-^F-^Wc)) =diYziF-^F-'^Wz{F-’^Wzc)). 

Proof. To obtain the first statement transform the well-known equation diVa;(cJ) = c. The 
second follows from A 2 ;(diva; j) = diVa;(A 2 ;j). Next transform diVa;(Curl 2 ;(c)) = 0 and div 2 ;(cj) = 
cdiv 2 ;(j) -I- Va: c • j. Finally observe that curl^-CVa; c) = 0 as well as curl 2 ;(A 2 ;j) = A 2 ;(curlj). 

If diva;(j) = 0, a simple calculation together with the fact that in this case = —§^ 
shows that Va;(curla; j) = [)* “g] Ax{j), which upon transformation yields the result. For the 
last statement consider S/^lA^c) = Ax{\/j;C). □ 

A.6 Remark. 

Let iy(x) be the unit normal vector at x € dll. Then the corresponding transformed unit normal 
vector is given by 

i/{x) = \\F~'^{x)v{x)\\~^ F~'^{x)v{x). 

If n = 2, the unit tangential vector f{x) has the direction ~q'\ F~’^{x)i'{x) = 

F{x) [? ~o] i^{x) = F{x)t{x), thus it holds 

•f(x) = ||A(a:)T(a;)|r^F(a;)T(a::). 

II'll indicates the chosen norm in K". 
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B Very weak solutions of the transformed Stokes system 

In this appendix we develop the theory of very weak solutions for the transformed Stokes equa¬ 
tions in Oi, which has been suggested for F = Id by Conca in |Con87| . We then derive estimates 
of the velocity and the pressure. We are looking for functions {B,j3) G x as 

solution of 


-div{F-^F-^B)+F-^V P = Gi+div{F-^G 2 ) in fli (B.la) 

div(F-^F) = 0 in fli (B.lb) 

B = ^ on Ey (B.lc) 

S,/3 are L-periodic in Xi. (B.ld) 


Here E^ := E U (0,F) x {h}, and ^ € L^(Et)^, Gi G as well as G 2 G are given 

functions. We require the compatibility condition 


- F ■ 62 da + 
JY. 


F • 62 dcr = 0. 


{X2=h} 


Define the space W 3 = {g G is L-periodic in Xi,J^^ g dx = 0} and let g G F^(r2i)^, 

u G W3. The auxiliary problem 


(B.2a) 

(B.2b) 

(B.2c) 

(B.2d) 


can be transformed to a Stokes system in Hi and solved on that domain. This yields the existence 
of a unique solution (</>, tt) G x with 


— div(F ^0) + F -K = g 

in Hi 

div(F“^0) = u 

in Hi 

0 

II 

on Yjx 

0 , TT are F-periodic in xi. 



II‘/’IL2(S2i) 2 + < <^(115111,2(12^)2 “k 11 U11 ^1 (j^^ , ) . (B.3) 

Using integration by parts, one calculates that 

j B ■ g dx — {(i,u)= j (Gi-k div(F“^G 2 ))<^ dx — j F~^{F~'^ \/(p — ttI)!/■ ^ da 

</fix J Yl'j' 

Define for {g,u) G x W 3 

l(g,u) = (Gi +div(F-iG 2 ), 0)121 - , 


where 0, tt is a solution of (B.2k 

B.l Definition. 

{B,j3) G X is called a very weak solution of the Stokes problem (B.l), if for 

all {g,u) G X IU 3 the identity 


{B,g)n^ - {(3 ,u)h~Gh^ = 


holds. 
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B.2 Lemma. 

The functional I : x W 3 


is linear and continuous. 


Proof. The solution operator to (B.2 1 is linear, as is the gradient and the duality pairing. Thus, 
one obtains the linearity of 1. For the continuity, we estimate 

\l{9,u)\ < |(Gi +div(F-iG 2 ),<^),J + 


L 2 (ni )2 II'/'IL 2 (J 2 j )2 + 11^211^2(5^^,4 11 ^^'11 „1 ( 5 ^ ^ ,2 + G || V ( 


< ||Gl 

^ (ll^llL2(ni)2 + 11^211^2(5,^,4) II<?!>L2(5,^,2 + ^ W^W 
^ C'(||Gi||,^2(5,^,2 + 11 G 2 11 ,,2 ( 5 ,^, 4 + 2 ||^||,,2(,,^, 

< C '(|| Gi ||,_ 2 ( 5,^,2 + || G 2 ||,_ 2 ( s ,^,4 + || C ||,^ 2(,,^)2 


ll, 2 (S 5 p)‘i ll^ll T C ||'^||,, 2 ,,,_^, 


lien 




lien 




■GUV' 


14 . 2 ( 5 , 1 ). 


lien 


I,2(S,,)2 


ff2(s2i)2 + 11 V '^111.2(5,1,) 


I, 2 (ni)- 


+ U 




o) 


by (B.3). This shows that I is bounded and thus continuous. 

B.3 Proposition. 


□ 


There exists a unique very weak solution of (B.ll. 

Proof. The Lemma above shows that I G (L^( 02 )^ x W 3 )*. Since L^{n 2 )^ x W 3 is a Hilbert 
space, the Riesz representation theorem yields the existence of a unique B G and a 

unique fi G IT 3 such that 

{B, 3 ) 5,1 + {P, u)ii-i_ 111 = l{g, u) V 3 , u G X W 3 . 


By choosing B = B, jS = —/3, we see that (H, /?) is a very weak solution of (B.l I. 
B.4 Lemma. 


□ 


Let {B,f3) be a very weak solution of (B.l). There exists a constant G > 0 such that the estimate 

(B.4) 


-^ 111 . 2 ( 5 , 1,2 < C ' dIG'l 111 . 2 ( 5 , 1,2 + 11 ^ 2111 . 2 ( 01)4 + 11 ^ 111,2 


(S.r)2 


holds. 

Proof. Choose g = B and u = 0. Using the estimates from the proof of the previous lemma 
together with the scaled Young’s inequality yields the result. □ 


C Various existence theorems 

C.l Transformed Stokes equation 


In this section we prove the existence and uniqueness of the solution of Problem (2.4) for fixed e. 
Basically, the approach is the same as in the functional-analytic treatment of the Stokes equation 
(see e.g. |SP80p . 

By multiplying (2.4aI with </> G iLjji,(U®)^ where 

:= {w G I div(F-iu;) = 0}, 

integrating by parts and noting that 

/ F~'^{x)\/p{x) ■ 4 >{x) dx = — i p(a:) div(F'“^(x)(()(a;)) da; = 0 , 
an® an® 
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we obtain the weak formulation of Problem (2.41 in the form 

f F~^{x) {x) ■. F~^{x) (l){x) dx = ( f {x) ■ (j){x) dx V(/) G (C.l) 

Note that is a Banach space with respect to the norm || V •||j^ 2 |.s^e) 2 . We need the 

following lemma for the estimation of the left hand side: 

C.l Lemma. 

There exist eonstants 0 < kp < Kp such that for the eigenvalues X{x) of F~^{x)F~'^{x) holds 

kp < X{x) < Kp Vx G fl, 

i.e. F~^{x)F~'^{x) is symmetric and positive definite. 

Proof. A calculation of the eigenvalues Ai(x),A 2 (x) of F~^{x)F~'^{x) yields 


Ai(x) — 1 + 






a'ixiY 


^ 1 , g'ixiY , g'{xiY 

X2{x) = 1 + — - - yg'Yir + — i —• 

Because of the smoothness of g there exists an M > 1 such that | 5 '(x)| < M for all x G fl. 
Obviously Ai(x) < 1 + 2M^ =: i = 1, 2. 

Choose a kp small enough such that + 2 < Another calculation shows that 


A2 (x) > kp 


g'{xiY < - \-kp-2, 


kp 


which gives the desired result since Ai(x) > A 2 (x). 

C.2 Proposition. 


□ 


Let e > 0 be fixed and let F be given by (2.3l. For given f G L^(0), the Problem (C.l) has a 
unique solution K G 

Proof. Define for u,v G the (bi-)linear forms 

a{u,v) = / A'“^(x) V u(x) : F“’^(x) V r!(x) dx 


and 


h{v) = I /(x) • v{x) dx. 


The continuity of 6 for / G Lf{LlY is standard. In order to apply the lemma of Lax-Milgram, we 
have to show that a is continuous and coercive. First note that as a pointwise estimate we have 

F~^{x) V v{x) : F~^{x) V v(x) = tr(V x(x)^F“^(x)F“^(x) V v(x)) 

2 

= ^ ef V x(x)^A“^(x)A“^(x) V v(x)ei 


i=l 

2 


< ^ ||Vx(x)e*||, ||F 1 (x)F '^(x)||j|Vx(x)ei| 


2=1 


< Kp'^2, ef V v{x)^ V v{x)ei 


i=l 


= KpVv(x) : Vv(x) = Kp ||Vx(x)||J 
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with II-11 2 being the Euclidean vector- and matrixnorm. This gives the continuity of a due to 


j |E ^{x) u{x) : F V u(x)| da; 

^ i j=i ^ 

= j (F“'^(x) V u(x) : F“^(x) V m(x)) ^ (F“’^(x) V u(x) : F“^(x) V u(x)) ^ dx 
JQ'^ 

<KFf ||Vu(x)||J|V?;(x)||, dx 


<ifFl|Vu|| 


IIV?;| 


where the Cauchy-Schwarz inequality in has been used in the last step. For the coercivity 
consider 


fcF||Vu||^2(„e) 


< / A 2 (x) V u(x) : V u(x) dx 

< / Xi{x)ef V v{x)'^ V v{x)ei dx 
Jne i=i 

< / F“^(x) V u(x) : F“^(x) V u(x) dx. 


Now the Lax-Milgram lemma implies the proposed result. 


Due to (C.ll, the solution it® fullfills 

- div(F-iF-^ V U®) - / e 


□ 


We will now characterize the orthogonal complement (i?j;y(f2®)^)-*~ of in order to 

reintroduce the pressure. We remind the reader of the following results, the proofs of which can 
be found in |Wlo92| : 

C.3 Theorem (Generalized Trace Theorem). 

Let k GM.,k > 2, and let K he a bounded domain in M”, n G N with boundary dA G . There 
exists a continuous linear operator T : H^{A) —> iJ^“2(9A) x F[^~^~"s [dA) with 

F{(j)) = {(j)\gA, ^IdA) for all (j) G (A). 

C.4 Theorem (Generalized Inverse Trace Theorem). 

Let A he a bounded domain in M", n G N, with boundary dA G with a given fc G K, fc > 2. 

Let T be defined as above. 

There exists a continuous linear extension operators : H^~"i {dA) x {dA) —> H^{A) 

such that 

r O £■ = Id. 
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C.5 Lemma. 

It holds 

Proof. Define G := {P-'^ \7p \ p G L^in^)}- Let cj) G G,u G with 4> = Vp. Then 

((ji,u)jj-i(fje )2 jji(fje )2 = (F Vp,u)jj -1 (fjs )2 jji(fjej 2 = / pdiv(F %i) da; = 0, 

an® 

such that (j) G Therefore G C 

For the other inclusion we will show that div(F“^-) : Hl{LP) —> L^{LP) is surjective and 
that —F“^ V • is its adjoint operator, therefore being injective from Fq(D®) to im(—F“^ V •). 
Now if '0 G we consider u G HQ(fl‘^) with div(F“^it) = 0. It holds 


(0,m) 




Since u is arbitrary, 

0 T ker(div(F“^-)), 

and since ker(div(F“^•))-*- = im(—F“^ V •) there exists a, p G F^(fl®) with 

tjj = F~^ V p. 


The surjectivity of div(F is a consequence of Lemma C.7 below; and the adjointness of 
the operators can easily be seen from the equation 


(F ^ Vp,'u)^-i(„e )2 ^i( 5 ,e ,2 = - / pdiv(F ^u) do; = (p, div(F ^■u))^ 2 (j 5 e,. 

an® 

Before proving some properties of the divergence operator, we need the following lemma. 

C.6 Lemma. 

Let 6 G {VP) ■ Then 

Curl(0) ■ V = — VJ 0 ■ T on dVP 
Curl(0) • T = V 9 ■ V on dVl^^. 


□ 


Proof. It holds 


Curl(6») • V = 


0 -1 

1 0 


v/e-u = v/0- 


(f? 


since the matrix [° corresponds to a rotation of | and thus [° o]^ ’ ^ second 


equality follows along the same lines. 

Now we are ready to prove the lemma used above: 

C.7 Lemma. 

Let G G L^(fl'^) with G = 0. There exists a (p G Hq{VP)‘^ with 

div(F“^(x)0(a;)) = G{x) in D® 
(j){x) =0 on dVl^ 


□ 


such that 


2 < 6 - 11011 ,, 




Thus div(F ^■) : — > LQ(r2^) is surjective. 


30 







Proof. We look for (f in the form 

(j) = FVr] + FCuTl{e) 

with 77 satisfying 

Arj = G in 
V 77 • 7 / = 0 on dOf. 

By considering the weak formulation of this problem 


f [ G ■ if 


and using estimates similar to those derived in Propositon |C.2| we see that a unique solution 
77 e i7p(0^)/R exists, satisfying the estimate ||V r]\\^ 2 ^^^s ^2 < C 
By regularity arguments one can show that 


9^77 


dxjdx.i 


<C\\F \\^2 






bJ € {1,2}. 


As for 6 , it should hold 


Curl(0) ■h' = —'\79-T = —'\/r]-iy = 0on 
Curl(6») -T = V0 - ly = -V Tj-T G on dn^. 


By the general inverse trace theorem C.4 there exists a 0 S H^{PF) with V 6 ■ lyldQe = — \7 rj ■ t 
and = 0 (thus especially V 0 • r = 0 on 90^) and 

||0|L.(,e,<C'||V77||,,^,.,. 

Now we have V 77 + Curl(0) = 0 on therefore also F{W 77 + Curl(0)) = 0 on the boundary of 




□ 


To reintroduce the pressure, notice that by equation ( C.l[ ) 
- div(F-iF-^ V 77^) - / € 


By Lemma C.5 there exists a pressure G Lf{FF), unique up to a constant, such that 

— div(P"“^(x)F“^(a:) V u‘^(x)) — f{x) = —F~'^{x) Vp‘^{x) 

holds in This finishes the considerations about the existence and uniqueness of the trans¬ 
formed Stokes equation. 

We have the following regularity result: 

C.8 Proposition. 

If f G r > 0, then G and G H''+^{n^). 

We do not give a proof, which can be carried out by adapting the regularity arguments for 
the usual Stokes equation (see e.g. |Tem77| '). For the interior of the domain, one can use the 
following argument: 

Applying div(F'“^-) to Equation (2.4al gives (by the second formula of Lemma A.51 

div(F-^E-^ Vp'^) = div(E- V) € 

Therefore p^ G where VF' is a strictly included subdomain of FF. Because of 

- div(E-iE-^ V 77 '^) = / - F-'^F/p^ G 
we conclude that G . 
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C.2 Boundary Layer Functions 

We define some function spaces that are used in the sequel: Let 

v = [z€ lL(^bl)" I V z e ^ e 

°° /n 


^ }, z is 1-periodic in xi| 


and 


Fdiv = {zG LL(^bl)' I V z e L\Zbi.)\ z e z = 0 on Q }, 

k=l ^ ' 

divy{F~^{x)z{y)) = 0, z is 1-periodic in xi|. 

Define W as the completion of Vdiv with respect to the norm ||z||j^ = ||V z\\^ 2 ^z ) 4 . The Poincare 
inequality in Z~ reads ||z||^ 2 (z -)2 < C||V z||^ 2 ( 2-,4 for all z £V. 


C.2.1 The Main Auxiliary Problem 

For the development of a theory for the boundary layer functions, we start with a more general 
formulation: 

Let 7 i > 0, CT G H^{SY, p G L?{ZY and pi G Lp’i^ZY be given. Assume that G 

L^(Fbl)^ and G L^(Zbl)^- Fix x £ fl and consider the following parameter-dependent 

problem: Find C G W such that 


f F ^{x)WyC{y) ■■ F ^{x)Wy(l){y) dy = [ p{y) ■ (l){y) 

'J Zbi. -^bl 


dy 


(C.2) 


- [ Pi{y)-F ^{x)^y(l){y) ^y+ [ F ^{x)a{y) ■ (l){y) day y <j) £ W 

C.9 Proposition. 

There exists a unique solution of Problem (C.2|. 

Proof. The result follows by application of the Lax-Milgram lemma: 

Define for u,(j) £ W 

B{u, (/))=[ F-^[x) Vy C(2/) : {x) (j){y) dy, 

K4')= [ p{y)-Hy)<iy- [ Pi{y) ■ F~'^{x)\/y(j){y) dy+ [ F-'^{x)a{y) ■ (j){y) day. 

-2'bL -^BL -S’ 

The continuity and coercivity of the bilinear form B inW can be proved analogously to the case 
of the transformed Stokes equation, see Proposition |C.2[ To see that b is bounded, note that 


p{y) ■ (fiy) dy 


< 


iz+ 


p{y) ■ (j){y) dy 


P{y) ■ </>(2/) dy 


< 11(1 + y2)p|li,2(2+,2 ||(i + 2 / 2 ) ^</>| 




\L^{Z-) 


-p\m 
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where we used the standard Poincare inequality in Z , the fact that |(1 + y 2 )p\ < < 

, see |,TM96] . The estimation 


I,2(Z+)2 — 


1 ^ 2(2 + ) 


□ 


^g7i|y2||p| < and ||(l + ?/ 2 ) ^(t>\ 

of the remaining terms is standard. 

C.IO Lemma. 

Let dYVy{F~^(x)pi(jj)) € L^(Zbl)^ o,nd let p,pi,a be 1-periodic in yi. Then the solution C, of 
( |C.2[ ) is in ■ 

C.ll Proposition. _ 

Under the assumptions of Lemma C.IO there exists a pressure field k G L^q^(Zbl) such that 

- diVy{F~^{x)F~'^{x) Vy C(2/)) + K{y) 

= p{y)+dWy{F~^{x)pi{y)) in W. 

Proof. We are going to nse analognes of Lemmas |C.5| and |C.7| for an increasing seqnence of sets 
in order to show that W-^ = {F~^{x) VyP \ p & -^?oc(-^bl)}- 

Define for I G N the sets Z* = [0,1] x ((0, 1) U (Ui=i{^* ~ (Di) space 

Wi = [zG H\Zff \z = diory 2 = ±l and on [j {dYg - Q }, 

z is 1-periodic in yi 


It is clear that Zf C Z*j^.^ and that each Zf is a Lipschitz domain. 

div;(F ~^(a: )-) : Wi — Lq{Z*), div/(F“^(a;)-) := divy(F“^(a;)-) is snrjective by an analogne 
thus F~'^{x)Vi{-) := F~'^{x)Vy{-) is injective from L^^Zf) to Wf 
€ V' snch that (/, ^Lzel)^ “ ^ ^ 


of Lemma 


C.7 


Now let / 

ker(div/(F“^(a;)-)) be given and denote by u the extension by 0 ontside Zf. Since then 
diVy(P“^(x)'u) = 0 in Zbl we have (/,'S)jf-i(z„L)^.~ duality of the extension 

operation we conclude that f\z* -L ker(div/(T'“^(a;)-)). Therefore f\z* G im(F~'^ (x) V i •), and 
there exists a pi G L^{Z*), unique up to a constant with / = F~'^{x) '^yPi in Z*. 

Since Zf G Zifj^.^, the difference pi+i — Pi is constant in Z* and we can choose p;+i in snch a 
way that pi+i = pi in Zf. Thns / = F“^(x) Vyp with p G ^^^^(Zbl)- 

The pressnre k can now be obtained by observing that - via an integration by parts of (C.2) 
dWy{F-^{x)F-'^{x)VyC,[y)) + p{y) + d:wy{F-^{x)pi{y)) G W-^. □ 

C.12 Lemma. _ 

Let ^ and k be defined as above. Under the assumptions of Lemma C.IO we have f G Hi^^{Z)^ 
and K G Hl^^{Z). 


Finally, we obtain the following strong form of Problem (C.2): 

- divj;(P"^(3:)F"^(a:) Wy ({y)) + P"^(a;) Vy /t(y) 

= p + diVy(F“^(a;)pi(j/)) a.e. in Z (C.3a) 

divy{F~^{x)C{y)) = 0 a.e. in Z (C.3b) 

C(2/i,±0) = on 5 (C.3c) 

C = 0on[jZAdys-O} (C.3d) 

C, K are 1-periodic in yi (C.3e) 


with known functions G H ^ (S")^. 
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C.2.2 Exponential Decay 

Define for k £ —N the sets Zk = Z~ n ([0,1] x [fc, /c +1]) (these domains, as well as other auxiliary 
sets needed in the course of the derivation, are depicted in Figure]^. 

C.13 Proposition. 

Let p := p + divy(F~^{x)pi) £ L^{Z~)^ and let C and k be as above. Define 

Then the following estimates hold: 




|rfc+i-r,|<C(|jVC|| 




+ 


li-Fzfc)' 


0 




0 


Proof. Define the space 

14 = jz G H^{Zuf I z = 0 on 9Zfc\(({0} U {!}) x [k,k + 1]), z is 1-periodic in ?/i|. 


Consider Equation (C.3al on Z^ with Vy(/t — r^) instead of Vy n. By multiplication with a test 
function and integration by parts we obtain 


f F ^{x)\'yCiy) ■. F ^{x)\'y(l){y) dy - [ {K-rk)diVy{F \x)(l){y)) dy 
J Z}^ J Z}^ 

= [ ip{y) + divy(E"\a;)pi(y)) • (/>(y) dy 'itf £ Vk. 

J Zv, 


Analogously to Lemma C.7 there exists (fk £ 14, solution of 

divy{F~'^{x)(l)kiy)) = K-Tk in Zk 


with 


I Vy <Pk\\, 


' V 4 C ||k ■ 


C depends only on the geometry of Y* but not on k. 

Inserting <f>k in the above equation and remarking that LF'“^(a;) Vy z 2 4 C |j Vy z ||^2 yields 


l«-^fclL=(z„ < 


'(’fell, 


|Vy<^fc|L 


II ■■ y A ^ ll^y ll-L^(Z(.)^ II '' y 'Xfelli2(2^)<i 

< C'(||p||^2(z^,2 + ||Vy CIL2 (z^,4) ||k - ^fclL2(z^, , 


thus the first assertion is proved. 

Next, set Zk^k+i = ZkU Zk+i and consider (fk,k+i satisfying 


divy(F 4a;)(()fe,fc+i(y)) = 


1 in ZO 

-1 in Z°^^ 


cfk.k+i = 0 on {dZk U aZfc+i)\(({0} U {!}) x[k,k + 2]) 

•Pk.k+i is 1-periodic in yi 

(the existence is assured since the right hand side of the first equation is in L‘^{Zk,k+i) and has 
mean value 0). 
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Figure 6 : Auxiliary domains based on the boundary layer cell Zbl- Left: The translated 
reference cells Z^, k £ —N. Middle: The two-cell subsets Zj^^k+i, illustrated with the sets 
Z-2,-1 (shaded with lines) and Z_3,_2 (shaded with dots). Right: The unbounded strips 
Z~{k). In the figure the sets Z~ {— 1 ) (shaded with dots) and Z~{— 2 ) (shaded with lines) 
are shown. 


Testing ( |C.3a[ ) with (l)k,k+i in Zk,k+i gives 

-[ K{y)dy+[ K{y)dy+[ {x)Vy(^{y) : {x)Vy (j)k,k+i{y) dy 

p(y) • (l>k,k+i{y) dy. 


/ Zk,k+i 


Note that ||(/)fc, fc-Hi 11^2,2^ 2 < C'||Vy(/i/,,fc+i||^ 2 , 2 _, ,4 < C\Zk,k+i\, thus dividing the equation 

by |y*| gives the estimate 


K+i-rfc|<C(||VCIL.<,,_,^^,4 

which finishes the proof. 




0 , 


□ 


C.14 Proposition. 

Por k G —N choose functions dk G C°°(M<o), 0 < dfe < 1 with dk{z) = 0 for z > k + 1 and 
dk{z) = I for z < k, z € M>o, such that dk and the derivative d'f. are bounded uniformly in k. 
For y = (^1) e [0,1] X (-oo,0] define akjy) := dk{y2)- 


Let ^,K be a solution of Problem (C.3l. Then it holds 


[ \P '^ix)\7yC{y)\'^akiy) dy = [ {K-rk)C,-P {x)Vy (Jk{y) dy 

Jz- JZk 

+ f Ky) ■ CiyWkiy) dy - [ p-^{x)VyC{y)p-'^{x){C{y) (^Vyakiy)) dy. 
J z~ J z- 
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Proof. Testing (C.3al with (f e C^(Zbl), </> = 0 on ~ (Di ^ 1-periodic in yi 

yields 

[ F-^{x)VyC{y) : F-'^{x)\/y(j){y) dy- [ Kdivy{F-^{x)(l){y)) dy [ p{y) ■ (j){y) dy. 

2'bL -2^BL -^BL 

Define for I <k — \ the functions ak.i = <Jk{l — (Ji). Choosing (p = fak.i leads to 


[ |T ^{x)WyC{y)\‘^ak,i dy = [ K{y)C{y) ■ F ^{x)Wy ak,i{y) dy 
Jz- Jz- 

+ [ p{v) ■ (.{y)^k,i{v) dy - [ F-'^{x)VyC,{y) ■. F-'^{x){C,{y)®Vy(Tk,i{,y)) dy, 

JZ- Jz- 

where we used the fact that 


F '^{x)Vy{C,(Jk,i) = F ^(a;)(Vy C)o'fc.i + C® T {x)y y <7k.i, 
diVy{F~^{x)Cak,i) = C, ■ F~'^{x)Vy Ok^i + CTfe,; diVy(T"^(x)C). 


We want to pass to the limit I — > —oo for fixed k. First observe that ak,i —> Ok as well as 
V ak^i — > V (Jk pointwise for I — > —oo. As |crfc ;| < C and | V ak,i\ = |(V a-k){l + <Ji) — ak V (7;| < 
C a.e. with a constant C, we obtain that almost everywhere 

I \F-^{x)Vyay)\^ak{y)\ < C\F-^{x)VyC{y)\^ 

\Ky) ■ Ciy)<^k,i{y)\ < C\p{y) ■ C{y)\ 

\F-^{x)\/vCiy) : F-^ {x){Ciy) (^kAy)) \ < C\F-^{x)VyCiy) : F-^{x){C{y)<g)I)\ 

where / denotes the identity matrix. Since the right hand sides are integrable, application of 
Lebesgue’s dominated convergence theorem yields for I —> —oo 


[ \F ^{x)VyC{y)\^(Tk,i dy 
Jz- 

p{y) ■ C{y)<xk,i{y) dy 


fZ- 


[ \F '^{x)WyC{y)\^ak dy 
Jz- 

/ p{y) ■ AyAk{y) dy 


and 


[ F ^{x)WyC{y) : F ^(x)(C(j/) O crfc,i(y)) dy 
Jz- 

—^ [ F-'^{x)\7yC{y) : F-'^{x){C{y)(g)Vyak{y)) dy. 

Jz- 

Finally we have to consider the term K{y)C{y) ' F~'^{x) Vy crk,i{y)- Because of V (Tk,i{y) = 0 
a.e. for y ^ Zk U Zi we have 

[ K{y)Ciy) ■ F~'^{x)\7yak,i{y) dy = [ (^(y) - rfc)C(y) • Vy crfc,i(y) dy 

J Z- J ZkUZ I 

= [ (Ay) - rk)C{y) ■ F~'^{x) Vy cjkAy) dy + / {Ay) - D)C(y) ■ f~'^{x) Vy (Jk,i{y) dy 

JZk Jzi 

+ {ri - Tk) / C(y) ■ F~'^{x) Vy (7k,i{y) dy. 

Jz, 
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For I —> —oo we obtain by using Poincare’s inequality 

(ri-rk) [ C,{y) ■ F~'^{x)\/y(Jk,i{y) dy 

Jz, 


<^^l|vcil 




and 


f -n)Ciy) ■ F '^{x)Vyak,i{y) dy 

JZi 


< c(||vcii 

^0, 




) l|vcil 




where we also used the preceding lemma for the last estimate. Thus arguing similarly with 
Lebesgue’s theorem one arrives at 

lim [ K{y)C{y) ■ F~'^{x)Vyak,iiy) dy = [ - rk)C{y) ■ F~'^{x)Vy ak{y) dy, 

1^-00 Jz- Jzu 


and the proof is complete. 
Define for k G —N the sets 


□ 


Z-{k) := Z” n ([0,1] X (-oo,fc]). 


C.15 Proposition. 

Let p G L'^{Z~)'^ and let (^,k he a solution of problem (C.3l. There exists a eonstant Cq inde¬ 
pendent of k sueh that 

11'^ ClL2(z-(fc), < Cq ||p|li2(2-(s,,)2 . 

Proof. We estimate the terms on the right hand side of the previous proposition separately: By 
the Poincare inequality 

[ F-'^{x) Vy C(y) : F-^(a;)(C(2/) 0 crk,i{y)) dy 
Jz- 

[ F-'^{x) ^yCiy) : F-^{x){C{y) 0 Wy ak,i{y)) dy 
JZu 


and 


<C||V.CIIz2,,^,4 

p{y) ■ C{y)crk{y) dy 


lz-{k) 


<c^l|v,CII 


2/ ^ II (fc))^ ll/^lll.2(Z“(fc))2 ’ 


p{y) ■ C{y)<^kiy) dy 
z- 

Using Proposition |C . 1 3] gives 

/ {k- rk)C ■ F~’^{x)Vyak{y) dy 

Because of Z~(k) C Z~(k + 1), Zk C Z~{k + 1) and Young’s inequality we obtain 

kp f I VyC( 2 /)Pcrfc(y) dy < / \F-'^{x)VyC{y)\'^(7k{y) dy 
Jz- Jz- 


<C'l|v,Cllz2<.,P + 




IVCII 




<C*\\^yal^^,,f + C5 [ |V,C(2/)Pdy+§ 

Jz-(k) ° 


\L^(Z-{k + 1))2 
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for (5 > 0. Next observe that 


[ I VyC(y)Pcrfe(y) dy = / I VyC(y)Pcrfc(y) dy + / \WyC{y)\'^(Jk{y) dy 
Jz- Jz-(k) JZk 


and 


[ |V,C(y)l"dy- / |V,C(y)l"dy 

Jz-{k+l) Jz-(k) 


thus leading to 

{kp -CS + C*) f I V, C(y)P dy < C* |1 V, Cll 

Jz-(k) 


Z-(k) 


L^(Z-(k + I))'* 


Cl 


\L^^Z-(k + 1))^ ■ 


Choosing S small enough such that kp — C6 + C* > 0 and kp > CS gives the recursion 

Ofe < jak+i + Fk, k € -N 


with 


Ofc = ||VC||^2,^-(,„ , 7 = 


C* 


kp-CS + C* 


< 1, 


Fk = 


Cl 


kp-C5 + C* + ■ 

Since Z~{k) C Z~{k + 1) we also have Fk < Cfe+i- This implies the claim as in |JM96| . 
C.16 Corollary. 

Consider the situation as above. Then there exists a constant Kao, 


□ 


Kao = lim 


1 


k—k—oo If* 


'Zk 


K{y) dy 


and a constant C*, independent of k, such that for k € —N holds 

k 


\K ^OO II j;,2|.2-{fc)) 


< C* ^ IIpII; 


L^(,Z-(l + 1))2 ■ 


1 — — 00 


Proof. Proposition |C.13| yields 


«:-rfc||,2,„,<C(||VC|| 


|rfc+i-rfc|<C(||VC|| 


r^iZkY 
r^(Zk,k + id ■ 


+ 




ip)- 


We show that is a Cauchy sequence in M, thus providing the existence of Koo- By the triangle 
inequality it holds for k G —N, I < 0 


\rk+i - r-fcl < ^ \rk+j - rk+j-i 
1 

3=1 

< C(|| V C||^2,2-(;, _ + ||plL2(2-(fc _ 1))2) 
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where the last constant is independent of k and 1. Since the last term converges to 0 for fc —)■ — oo, 
we obtain the desired resnlt. 

Next observe that 


1^ ^oo||r,2,z„i ^ ||k II + Ir m ^oo|li,2(z„) 


— II '"Hi' (Zm) 

<c(||vciL2(,^,4 + ||pIL.(,„, 

< C'djV CIL2(2-(„+ 1))4 + ||pIL2(z-(„ + i„2) 


) + |y*| {r^ + Kooi 

(Z-(™ + 1))2 


+ |y*| lim 

j—^oo 

E ’’m-i “ ^m-(i+l) + U-(i + l) 
^ 1=0 

^00 1^ 

<c^IIpIL7z-o 

. + .„. + i:c(l|vciL., 

l-o 

,4 + IIpII s^2) 

<c^IIpIL7z-o 

^ + 1))2 + 2 C'(||VC||^ 2 ( 2 -(„ + i,)4 + 

ll^lli 2 (Z“{m + 1 ))^) 

<c^IIpIL7z-o 

T. + 1))2 I 


the above inequalites and Proposition C.15 
— «^oo| = 0. Thus by 

Furthermore, note that 


1^ ^oo II j;,2(2-(fc)) 




m——oo 


<c- y: 


ll,2(z-(„+ i))2 


m— — oo 


the second assertion holds. □ 

Finally we are able to get a resnlt on the decay of the solutions C, n in the porous part Z~ of 
^bl: 


C.17 Corollary. 

Assume that p G L^(Zbl)^ for a 71 > 0. 
solution C, K of Problem (C.3) holds 


Then there exists a fi > 0 such that for the 


11^ CIL2(2-(fc„4 

II Cll i,2(z- (fc))2 

||k — Koo|L2(2-(fc)) 




Proof. By the assumption on p note that ||p||j^ 2 ( 2^)2 < Ce 'I'dfel, Therefore 


0 


L2(,-(,„2<C ^ ^ 

k— — oo k——(X) 


1 - e'y^ 


where we used the formula for the geometric series for e^^ >1. Using the same argument once 
again, one obtains 

E llplL2(,-,„<Ce-^^l'=l, 

1 — — OC 
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which gives the first and the last assertion. The second one follows due to Poincare’s inequality. 

□ 

In order to deal with the behavior of C and k in Z~^, we are going to use the theory for 
the exponential decay of solutions of elliptic problems, developed by Landis/Panasenko and 
Oleinik/Iosif’jan, see |LP85| and |OI81| . 

C.18 Theorem (Exponential Decay). 

Let the geometry he given as above. In Z'^ consider the elliptic equation 


-divy{F{y) Vyu(y)) = f{y) 

with a given matrix function F S L°°{Z'^)‘^ satisfying the following ellipticity condition: Let 
there exist constants ci^Ci > 0 and M > 0 (M = 1 in case F is symmetric) such that for all 

ci\e <eF{y)f<C^\e 

\v^F{y)^\ < M{if F{y)r])^ F{y)£,)^. 

Assume further periodic boundary conditions on ({0}U{1} xK>o) and Dirichlet and/or Neumann 
conditions on S such that there exists a solution u with V it G L‘^[Z'^). Let there exist constants 
q,Q> 0 such that f G 

Then there exist constants qi,Qi > 0 and Cu such that 

11^ ~ C'u||^ 2(2+ n ([0, 1] X [fe, oo))) — ■ 

Furthermore, there exists y* > 0 with 


\u{y) -Cu\< Qie for y2> y*- 

Proof. Theorem 10 of [OI81| gives the first two estimates. 

Due to the lifting property of elliptic operators, we obtain a solution u € H^{Z~^)-, and because 
of the embedding H‘^{Z~^) ^ C^(Z~^) there exists a continuous representative. Therefore we can 
apply Theorem 2 in |LP85| in order to get the pointwise estimate. □ 

C.19 Proposition. 

Assume that (/ k) is a solution of Problem 
and diVy{F~^{x)pi) G H^{Z^/. 

There exist / > 0, y* > 0, a vector Cq G and a constant such that 


C.3I with G H^{Z+f, e'^^y^pi G H\Z+/ 


11^ Cl 1 1 , 2 ^ 2 + n ([ 0 , 1 ] X [fc, c 

»))p - 

lie ^cIId2^2+ n ([0, 1] X [fc, c 

-)))^ - 

II ^ II n ([ 0 , 1 ] X [fc, c 

)) - 


Kiy) - c/< ce-^y^ 

for 1/2 > y* 

Ky) - C/< Ce-^y^ 

for 1/2 > y* 
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Proof. Set f = curl(C). By taking the curl of Equation we obtain due to Lemma |A.5| 

diVy{F~^{x)F~'^{x)\/y^{y)) = - cuTl(^p{y) + diVy{F~^{x)pi{y))'^ in Z+. 

The right hand side decays exponentially, thus by the preceding theorem we obtain 


V (curie) 

curl e — Cc 


i2(Z+ n {[0, 1] X [k. oo)))2 


< Ce-^^ 


L^(Z+ n {[0, 1] X [k, exj))) 


< Ce 


-pk 


with some constants j3 > d and Cc- 

Using the 7th assertion of the transformation lemma [ATS] we see that 


F-^{x) Vy{i{y)) = F-^{x) V,(curlC( 2 /)) 


0 

1 


-1 

0 


divy{F \x)F ^(x) VyC(2/))- 


Therefore 


divy{F \x)F ^{x)VyC{y)) =h{y) in Z+, 

h being a known function with € L^(Z~^). Theorem C.18 now shows that the first two 

asserted inequalities about the decay of ( and V ( hold. 

By taking the transformed divergence of Equation ( |C.3a[ ) one obtains 

divy(F~^(x)F~'^(x)VyK(y)) = - diVy(^F~\x)[p(y) + divy(F~\x)pi(y))]J in 

Again, the right hand side decays exponentially, and the estimate for k is proved. The remaining 
two inequalities follow easily. □ 

At the end of this section, we want to obtain some information about the constant Cr'. 


C.20 Lemma. 

For the solution ^ of Problem (C.3) it holds for all z < 0 


f ({yi,z)-F '^(x)e 2 dyi = f ({yi,0) ■ F ^(x)e 2 dyi = 0, 
Jo Jo 


and JgC ■ F ^{x)e 2 day = 0. 

Proof. By density it is enough to show the claim for () G VF n (Zbl)^. 

Integration of the equation divy{F~^{x)f{y)) =0 over [0,1] x (z, 0) and application of Stokes 
theorem yields due to the periodic boundary conditions 


f C{yi,0)-F ^(x)e 2 dj/i = f ({yi,z)-F {x)e 2 dyi =: Vz < 0. 

Jo Jo 


Since C ^ W it holds 



(^J C(yiU) ■ F ^(x)e2 dyi^ 


dt < C 




IC(yiU)P dyi dt < oo, 


and thus = 0. 


□ 
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C.21 Lemma. 

For the constant it holds C(^ ■ F~'^{x)e 2 = 0. 

Proof. Arguing as in the proof of the above lemma, we obtain 


[ C{yi,k)-F ^(a;)e 2 dj /1 = [ C{yi,-k) ■ F ^(a;)e 2 dyi 

JQ JO 


for all fc > 0. Now the left hand side converges exponentially to Cq ■ F '^{x)e 2 dyi = 
C(^ ■ F~"^{x)e 2 , whereas the right hand side converges to 0. □ 


C.2.3 Application to the Stokes Boundary Layer Problems 

We apply the results of the foregoing section to the problem: Find G V x Lj^q^(2'bl) 


such that for fixed x G it holds 

-divy{F-^{x)F-'^{x)Vyl3'°\y))+F-'^{x)\/yUj'°\y)=0 inZ (C.4a) 

divy(p-i(a;)/3“(z/)) =0 in Z (C.4b) 

[/3“]s(y) =0 on S' (C.4c) 

[{F-\x)F-^{x) Vy /3“ - F-\x)uj'°^)e 2 ]s{y) = on S (C.4d) 

^“(y)=0 on[jZAdYs-{l)} (C.4e) 

w'’* are 1-periodic in j/i (C.4f) 

where K^^{x) = F~^ {x)F~'^{x)V vP{x)e 2 . This corresponds to the case p = pi = 0 and a = 
K'°^{x). Lemma C.21 shows that 

C“-P"^(a;)e2 = 0. (C.5) 

Finally, we can obtain the complete information about the constants: 

C.22 Lemma. 

For all 0 < a < b it holds 


f w“(2/i,a) dyi = [ uj'^\yi,b) dyi. 

Jo Jo 

Thus the constant arising in the stabilization of the pressure is given by 

= [ w“(yi,+0) dyi. 

Jo 


Proof. Due to (C.4bl and the actual entries of F it holds 


(C.6) 


(C.7) 


Note that F ^{x) Vy = divy(F ^(x)a;°^(y)) (cf. Lemma A.51, thus Equation (C.4al reads 


column-wise 


diVy^F ^ix){F Vy/3S’'(y) - a;“(y)ei)) = 0 
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and 


diVyi^F \x){F ^(a;)Vj^/3^'(y)-w“(y)e2)) =0. 

Let 0 < a < b. Now integration of the eqnation 

0 = divy(j^-i(a;)(F-'^(x) \7y /3^'(?/) - a;“(?/)e2)) 

- g'{xi) divy(^F-^{x){F-'^{x) \/yP'^\y) - a;“(?/)ei)) 


over the rectangle [0,1] x [a, b] yields due to Stokes’ theorem and the periodicity of and 
in j/i-direction 

0= / (F"^(a;) - w“e 2 )(yi,&) • F"'^(x)e 2 dyi 

Jo 

1 

{F~'^{x) Vy / 3 “ - w“e 2 ) (yi, a) • F~'^{x)e 2 dyi 

-g'ixx) [ (F"^(a;) Vy/3i'-a;“ei)(yi,5) • F"^(a:)e 2 dyi 
Jo 

+ 9 '{xi) [ (F"'^(a;) Vy/3i'-a;“ei)(yi,a) • F"'^(a:)e 2 dyi. 

Jo 

Now we have (by using the actual form of F~'^{x)) 

w“(y)e 2 • F~'^{x)e 2 = w“(y) and a;“(y)ei • F~'^{x)e 2 = -y'(a;i)a;“(y) 



as well as 


F {x)V y PWy) ■ F '^{x)e 2 = {l+g\xxf)-^^p\\y)-g'{xi)-^^l3\\y) 
F-^{x) \/y /3'^\y) ■ F-'^{x)e 2 = (1 + g' {xif)^^0f{y) - g' {xi)^^0f{y) 




d 


dy2 


dyi 




where the identity (C.7) was used in the last equation. Substituting these results in the above 
equation, one obtains 

0 = ^ (^-g'{xi)-^P^^ - {1 +g\xif)-^P^^ + g\xi)-^P^^yyi,b) dyi 

-{l+g'{xif)f a;“(yi,&) dyi + (1+ y'(a;i)^) / u;“(yi,a)dyi. 

Jo Jo 

The first two integrals vanish due to the fundamental theorem of calculus and the periodic 
boundary conditions. We divide by (1 + g'{xi)^) to obtain 

f w“(yi,a)dyi= [ w“(yi,6)dyi V 0 < a < 6. 

^0 ^0 
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This proves the first statement. 

To obtain the second one, notice that for k > 0 due to Jensen’s inequality 


w (?/i, 2/2) dyi - C, 


bl 


dy2 < / |a;“(?/i,?/2) - dy 

J [0,1] X (fc,oo) 

—> 0 for fe —>■ oo 

because of the exponential stabilization of therefore dyi converges to C^* for 

6 —>■ oo. Now letting a —>■ +0 yields the result. □ 

C.23 Lemma. 

For the constant C'^* appearing in the exponential stabilization of the velocity it holds 

■ F{x)ei = f (3^\yi,F0) • F{x)ei dyi. 

Jo 

Proof. Let 6 > 0. Similarly to the above lemma, we multiply the equation 
0 = diYy(^F~^{x){F~'^{x)Vy l3i\y) - w“ei)^ 

+ g'{xi) divy(^F-^{x){F-'^{x) Vy (dfiy) - a;“e 2 )) 
by 2/2 and integrate over [0,1] x [0, b]. Integration by parts then yields 


(C.8) 


0 = - 


' [0,1] X [0.6] 


F ^{x){F Vy/ 3 “(y)-a;“(y)ei) • 62 dy 


-g'ixi) f F ^{x){F '^(x) Vy/SaHy) - w“(y)e2) • 62 dy 

a [0.1] X [0,6] 

+ f Vj;/ 3 “ - u;“6i)(yi,6) • F“^(x) 62 dyi 

Jo 

+ y'(xi) f b{F-^{x) Wy - co“ 62) (yi, 6) • T’-^(x)62 dyi. 

Jo 

As in the proof of the preceeding lemma we have T'“^(x)u;'^^(y)62 • 62 = a;'’^(y)62 • F“^(x)62 = 
a;’^'(y) and w*’'(y)6i • F~'^{x)e2 = —g'{xi)uj°^{y), thus the terms containing the pressure 
cancel out and we have 


0 = - 


' [0,1] X [0.6] 


F ^[x)F (x) Vy/ 3 f‘(y) • 62 dy 


-y'(xi) [ F '^{x)F ^(x) Vy/ 3 ^^(y) • 62 dy 

J[0.1]x[0.6] 

+ [ 6F"^(x) Vy/ 3 i'(yi, 6 ) • F"^(x)e2 dyi 

+ 9'{xi) [ bF~'^{x)\'y P2'iyi,b) ■ F~^{x)e2 dyi. 
Jo 
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Another integration by parts of the volume terms now yields 

0 = - / /3i\yi,b)F~^{x)F~'^{x)e 2 ■ 62 dyi 
Jo 

+ [ / 3 Y{yi,+ 0 )F~'^{x)F~'^{x)e 2-62 dyi 

-g'ixi) [ l32\yi,b)F~^{x)F~'^{x)e2 ■ 62 dyi 
Jo 

+ g\xi) [ /3“(?/i,+0)F"^(a;)F"^(a::)e2 • 62 dyi 

+ [ bF-'^{x) Vy I5\\yi,b) ■ F"^(a:)e2 dyi 
Jo 

- g'{xi) [ bF-'^{x) Vy I32\yi,b) ■ J^"^(a;)e 2 dyi. 

Jo 

When passing to the limit b —>■ 00 , the last two integrals vanish since Vj, decays exponentially 
to 0. The terms and converge to and repectively. Thus 

J (Ci' + 5 '(a;i)C' 2 ') A"^(a;)F"^(a;)e2 • 62 dyi 

= J (/^i' + g'{xi)f3^^'j {yi,+Q)F-^{x)F-'^{x)e 2 ■ 62 dyi. 

Since F~^{x)F~'^{x )62 -62 = l + y'(xi)^ we can divide the above equation by l + y'(xi)^, leading 


to 


CfFg'ixOCf = /'(/??' + y'(xi)/ 32 “)(yi,+ 0 ) dyi. 
Jo 


This is equation (C. 8 I. 


□ 


C.2.4 Dependence on the parameter x 

We summarize the results about the decay of the boundary layer function in the following propo¬ 
sition. Here, we take the dependence on the parameter x explicitly into account. 

C.24 Proposition. 

For the decay function it holds C^'{x) = fi^^x, yi, -fO) dyi. 


Proof. With the help of ( |C.5 1, (C. 8 I and Lemma C.20 it is possible to obtain the value of C'°^\ 
Using the exact form of F[x) and F~'^{x), the above conditions read 

Cf + g'{xi)Cf = PWyi.+d) dyi+g'{xi) P^\yi,+Q) dyi 


-g'{xi)Cf + Cf = -g'(xi)fg /3S’'(yi,-fO) dyi + /3^'(yi,-f0) dyi = 0. 


Since the determinant of the coefficient matrix fulfills 

1 g'(xi) 
-g'(xi) 1 


det 


= (l+y'(xir)^0, 


we can multiply both sides from the left by 


1 g'(xi) 
-g'(xi) 1 


-1 


to obtain the result. 


□ 
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By using the implicit function theorem for Banach spaces, one can show that the boundary 
layer function is continuously differentiable in the parameter x. By the Proposition above, 
this also leads to C'°'" having the same property. Denote by and the derivatives 
and for z = 1,2. They fulfill the equations 


-divy(F ^{x)F y) 

= d\Vy{-^^[F~^{x)F~'^{x)]\Jyl3'°\x,y))-^^F~'^{x)\IyUj°\x,y) m. Z (C.9a) 

divy(i^"^(x)^“’^(x, 2 /)) =-divj^(^F"^(a;)/3“(x,y)) in Z (C.9b) 

[/3“’®(a;)]s(y) = 0 on S' (C.9c) 

l(F-^(x)F-^(x) Vy - F-\x)u;^^’\x))e2]s(y) = 

-l^lF~^(x)F~'^(x)]Vy/3'^\x)-^F~^(x)uj'^\x)]s(y)e2 on S (C.9d) 

OXi uXi 

/3“’i(x,y) = 0 on - (Dl (C.9e) 

/jbi,i^^bi,i 1-periodic in yi (C.9f) 

as well as 

— divy(F~^(x)F~^(x) Vy 0°^''^{x, y)) + F~'^{x) Vy y) = 0 in Z (C.lOa) 

divy(F“^(x)/3'^^’^(x, y)) = 0 in Z (C.lOb) 

[/3“’^(a;)]s(y) = 0 on S (C.lOc) 

[[F-\x)F-^{x) Vy _ F-\x)J^^’\x))e2]s{y) = on S (C.lOd) 

/3“’2(a;,z/) = 0 on - (^1 (C.lOe) 

/jbi, 2 ,^bi ,2 1-periodic in yi (C.lOf) 


dne to i(a;) = 0. Since gf^F ^(a;) 


-9''{xi) 

0 


we have that 



A 

'dxi 


divy( —-F" ix)/3 ) 



d 

dyi 


A dy = 0 


due to the periodic boundary conditions for Thus a variant of Proposition C.25 below shows 
that there exists a function 0^ € Fl^^^^Z^n), Vy 9^ G F^(Z'bl) such that 


diyy{F-\x)9P{y)) = -diVy( —F-i(x)/3“(x,y)) 

in Z 

(C.ll) 

e^{y) = 0 

on Ur=i{^Vs - (°)} 

(C.12) 

[As(2/)=0 

on S 

(C.13) 

9^ is 1-periodic in yi. 


(C.14) 


Using 9^ to correct the divergence in equation (C.9b), one can use the results for the boundary 
layer functions to obtain an exponential decay towards constants C''’*’® (for and (for 
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w*^*’*), and the identities 

C'“’*(2;)=y /3“’*(a;,2/i,+0) dyi = 0°\x,yi,+Q) dyi =-^C'^\x) 

= y w“’*(a;,yi,+0) dyi = -^ uj'°\x,yi,+0) dyi = 

hold for i = 1,2. This shows that the derivative in a;-direction of P^^{x), de¬ 

cays to the corresponding derivative of the decay function C^^{x), C^^{x), and terms like 
S/xif^^^ix.y) — H{y 2 )C^^{x)) show the same decay behavior as li^^{x,y) — H{y 2 )C^^{x), lead¬ 
ing to similar estimates. 

C.3 Functions for the Correction of the Divergence 

In this section we consider the auxiliary problems associated with the correction of the trans¬ 
formed divergence Fix a; G fl and define 

Cp{x) = F{x)(^ f diVx(^F-\x)[/3'^\x,y) - H{y2)C'^\x)]^ dy^ea. 

C.25 Proposition. 

The problem: Find 9 such that 


diVy{F-\x)e{y)) = dWx(^F-\x)[^'^\x,y) - H{y2)C'^\x)] 

e{y) = 0 
[9]s{y) = C^{x) 

9 is 1-periodic in yi 


in Z 

on ur=i{^^s- (Dl 

on S 


has at least one solution 9 G F[^{Z)^ C\C^^[Z). 

Proof. We argue similarly to Lemma |C.7| and carry out the following ansatz: 

9{y) = F{x) Wy y{y) + F{x) Curly ^(y), 

where for y it holds 


in Z 


^yViv) = divx(F~^{x)[l3'^\x,y) - Ft(y2)C'“ 

Vyy(y)-!^ = 0 on U^i{9Fs - (°)} 
[Vy y(y) • e 2 ]s(y) = /" div,r(-F"^(a::)[/3“(x, y) - iL(y 2 )C'“(x)]) dy on ^ 

Ms(y) =0 on s 

y is 1-periodic in yi. 

We investigate solvability in the space Wd /K, with 

Wd = {z€ L'^^^iZBh) I Vz G L^(Zbl), -z is 1-periodic in yi}. 
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Define the linear functional 


C{ 4 >)= / div^(J^ ^{x)[( 3 '°\x,y)-H{y2)C'^\x)]) <i){y) dy 

-/ div:j,(^"^(a;)[/?“(a;,?/) - i7(?/2)C“(a;)]) d2/)(/)(yi,0) dyi. 


Since £( 1 ) = 0 the linear functional is well defined on Wij/K, and by the properties of y) — 

H{y2)C^^{x) it is continuous. An integration by parts shows that the weak formulation of the 
above equation reads 


■^BL 


Vy 77 • Vy (/) dy = £{ 4 >). 


Thus we get a solution 77, unique up to a constant. 
Next, we search for ^ satisfying 


d£ 

Curl (0 • 7/=-^ = 0 

oyi 

Curl(^) • T = ^ = - V 77 • T 
OV 2 


on ur=i{a>^s 
on ur=i{a>^s 


©}• 


Application of the inverse trace theorem C .4 to each cell Y — (°) in Z 


yields the existence of ^ similar to the proof of Lemma |C. 7 | 


and setting ^ = 0 in 

□ 


Since the right hand side of the equation for 77 decays exponentially, we can apply Theorem 
|C. 18 | and obtain an exponential stabilization of 77 towards some constant and a stabilization of 
V 77 towards 0 . As the construction of ^ is local, the decay carries over to this auxiliary function 
as well, and we obtain an exponential stabilization of 0 to 0 in 7/ for I7/2I —> 00. 

Therefore we obtain 


C.26 Proposition. 

The above problem has at least one solution 6 G V such that there exists a 70 > 0 with 

g7o|y2ly g 
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